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UNSTEADY FLOW IN CONDUITS WITH SIMPLE SURGE TANKS 


eH. Taylor,1 A.M. ASCE, Arnold Reisman? and Jack W. Ward? 


ABSTRACT 


The basic principles of the operation of simple surge tanks have been 
known for decades. American and European literature contains many articles 
devoted to this subject. In spité of the fact that the fundamental flow equa- 
tions are simple in concept, they are non-linear and hence difficult of solu- 
tion. The bulk of the existing literature with one relatively recent exception(1) 
has been concerned with numerical and graphical methods of solution. This 
paper presents a general solution of the simple tank problem in terms of a 
Single dimensionless parameter for the case of instantaneous and complete 
load rejection. The mechanical differential analyzer of the University of 
California was used in obtaining the solution. Seventeen curves covering the 
range from the case of practically no damping to that of almost critical damp- 
ing are presented. Whereas nothing new in principle is presented, it is be- 
lieved that the remarks below constitute a refinement of existing solutions 
and present numerical data covering a wide range of conditions in concise 
manner. 


INTRODUCTION 


The problem under consideration involves the prediction of the water sur- 
face elevation in the surge tank, depicted schematically in Fig. 1, as a func- 
tion of time following rapid and complete cessation of flow in the conduit. 


Note: Discussion open until July 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1933 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 2, February, 1959. 
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HY 2 UNSTEADY FLOW 3 
Theory 


The fundamental equation of unsteady flow in the conduit, based on the so- 
called “rigid water column” analysis, is, in consideration of Bigs: 


LW (1) 
H-h~ Heriction = g dt 
2 
- L@ ellos Vie (2) 
or z = eS ems 


in which z is measured positively upward. Very simply interpreted, the equa - 
tion states that the driving head less the friction head is equal to the inertia 
of the water column in the conduit. It is assumed that the inertia and friction 
of the water column in the tank are negligible in comparison with those in the 
pipe. 

The continuity condition at the conduit-tank junction requires that 


av ce bi (3) 


for the case of instantaneous and complete closure. 
Differentiating (3): 


dv A d“z 
Sige Sane (4) 
V and wv are eliminated between (3), (4) and (2) upon arrangement of terms: 
2 2 
an ad vs ae ss ee 
OF fetes eae ae ral Re = 0 
simplifying: 


d@z , f A (de\ ga, _ 
fy +84 +e. : (5) 


To account for damping during flow reversal in the conduit it is necessary to 
write eq. (5) in the following form. 


d@z , £ Adz|dz/ , ga ; 
Hts lal +h = 0 (6) 


This merely precludes the possibility of the damping term having opposite 
sign from the other terms in which case energy would appear to be added to 
the system by dissipative forces. 

Equation (6) may be rendered dimensionless by the following changes of 
variables: 
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Let 


T 

where 

tank system. 
Further let 


It is worth noting that T is the period of undamped oscillations in the conduit- 


¢4 
Zo 
in which Zp) is the driving head for steady flow 


d 
Then atu a Zo at 
a wi BE gcd Vabinan LAL. 
a a ~ do & dt 
dz ast 
(a) a = °° 
2 2° 
2s a2¥ ag 1 
b es 
(b) ae = 40 so at = 
Equation (6) now becomes 
Oo ecc din heth cad-sih AOL 
Zz ae 4 +tAs GI) S91, gae 
° a0 e485" hala kee 
Multiplying by %o°/2o 
a25 | fas, ap 
at * Ee 


a 
da do a3: (2m)? 5 = 


0 
The dimensionless coefficient 


AZo 
2da ’ 
the damping factor for the normalized equation. Finally 


ay 4 
aoe Ve 


further designated by the symbol wy , is 
ae + (2m)2? = 0. 


The initial conditions, corresponding to the instant the valve is closed are 


(8). 
_ ot a2) 
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Computer Solutions 


A diagram of the differential analyzer mechanization is shown in Fig. 2. 
Our remarks will be restricted to this particular problem since the theory of 
jthe differential analyzer has been presented in the literature.(2,3) For the 
most part, standard techniques were used in obtaining solutions to Eq. (8). 
The equation was written in the form 


a2 5 ea pas] ery2 
a = - ps [Se - (277) y 


and the right-hand side of this equation constructed by repeated integration, 
constant multiplication, and summation. The use of an integrator (No. 2) to 
multiply by a constant is standard where a large number of parameter values 
are used. A clutch or “disconnect” between integrators 3 and 1 in order to 


introduce a factor of -1 whenever af became negative (so as to obtain 


) was the only special feature involved. Seventeen solutions to eq. (8) 
as determined by the above technique are presented in Figs. 3,4 and 5. A 
dimensionless plot indicating the height and time of the first maximum surge 


is given in Fig. 6. 


CONCLUDING REMARKS 


It was found possible to express the differential equation describing the un- 
Steady flow in a surge tank system during a sudden outflow cessation in terms 
of dimensionless variables and a single dimensionless “damping” parameter. 
This equation was solved mechanically and the results are presented in a 
form which makes possible the determination of the maximum surge, the 
time of its occurence and an approximation to the time of complete cessation 
of water pendulation, for most simple surge tank system configurations. 

The computer solutions were substantiated within 1% accuracy by numeri- 
Cal analysis of a typical system.(4) An experimental model study of the 
University of California, Los Angeles (unpublished) likewise indicated sub- 
stantiation of the computer solutions to the extent that the height and time of 
the first maximum surge were found to be in excellent agreement with the 
above theory. However, the rate of damping as observed in actual flow has 
been tentatively found to be at variance with that predicted by the theory. In 
addition, the conditions of stability of a simple tank system were found ex- 
perimentally also to be in disagreement with the Thoma criterion which is 
widely quoted in the literature.(5,6) These latter matters are currently sub- 


ject to investigation. 
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EVAPORATION OF LAKE ONTARIO 2 


Ira A. Hunt, Jr.1 


SYNOPSIS 


There is definite evidence that the precipitation falling over the surface of 
the Great Lakes is considerably less than the precipitation falling on adjacent 
land areas. Evaporation computed by the water budget method is reduced in 
the same amount as the precipitation. The annual amount of net evaporation 
from Lake Ontario computed by the water budget and mass transfer methods 
is of the order of 2.0 feet. 


INTRODUCTION 


The advent of the St. Lawrence Seaway projects has stimulated the diverse 
interests of landowners, shippers, economists, and engineers, all of whom 
have been concerned about the levels of the Great Lakes for some time. Yet, 
despite their attentive concern, many are not aware that the outflow of the St. 
Lawrence River and, thus, the water levels of Lake Ontario, have been com- 
pletely regulated by the Iroquois Dam control structure since 1958. 

It is the responsibility of the Corps of Engineers, U. S. Army, to ascertain 
for the United States that the regulation of Lake Ontario provides optimum 
benefits for navigation, power, and riparian interests. To this end, the U. S. 
Lake Survey has been studying techniques of forecasting supplies to the Great 
Lakes for future use in effecting the most beneficial regulation plans. These 
studies determined that, of all the hydrologic factors affecting the levels of 
the Great Lakes, evaporation was the factor about which the least was known. 

Evaporation is the process by which water is converted to vapor. The 
molecules that make up water are in continuous motion. Some molecules of 
water have sufficient momentum to break through the surface and enter the 


Note: Discussion open until July 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1934 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 

_ American Society of Civil Engineers, Vol. 85, No. HY 2, February, 1959. 


a. Presented to the February 1958, ASCE Convention in Chicago, Ill. 
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air. The rate at which these water particles leave the water and enter the 
adjacent air depends upon the heat supply of the water and the condition of the 
air. Concurrently, molecules of water vapor in the air strike and enter the 
water; that is, they condense. In hydrology, evaporation is considered to be 
the net loss of water mass. 

There are several methods of computing evaporation. For the Great 
Lakes, the water budget method has been utilized most frequently. More 
recently, evaporation has been correlated with the vertical distribution of 
moisture in the air and the intensity of turbulent mixing. This is the so- 
called mass transfer theory. Another, called the energy budget method, 
takes into consideration solar radiation, atmospheric radiation, advected 
energy, sensible heat transfer, and other factors which affect the total 
amount of heat available to convert water into vapor. Computing evaporation 
by the water budget and mass transfer methods will be discussed herein. 


Water Budget Method 


The water budget method of determining evaporation on Lake Ontario is 
simply a matter of solving the general hydrographic equation for evaporation: 


Ei go yl PR, tt Fy ig HiT oe Oy gt By (1) 


where 


E_9 is the evaporation from the surface of Lake Ontario during the inter 
val 1 to 2. 


S,-2 is the change in the mean level of Lake Ontario during the interval 
1 to 2. 


Ry _9 is the runoff of water from the land surface draining into Lake 
Ontario during the interval 1 to 2. 


Py 9 is the precipitation falling directly on the surface of Lake Ontario 
during the interval 1 to 2. 


I,_9 is the inflow to Lake Ontario from Lake Erie, or the flow of the 
Niagara River plus the Welland Canal, in the interval 1 to 2. 


01-2 is the outflow from Lake Ontario, or the flow of the St. Lawrence 
River, in the interval 1 to 2. 


U,j_9 is the sum of all factors affecting the rise or fall of the level of Lak 
Ontario not contained in the other factors of the equation. 


Although this method of determining evaporation is simple, the determina. 
tion of the component factors which make up the equation is quite difficult. I 
studying Lake Ontario, the amount of hydrological and meteorological data 
available in the period prior to 1934 was limited, so a twenty-year period of 
record from 1934 through 1953 was used. Evaporation was computed in 3 
monthly units because the majority of data was available on a monthly basis 
and because it was intended that the results of the study would be used in 


regulation plans for Lake Ontario where a thirty-day forecast of the water 
levels is desired. 
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} Normal changes in levels of Lake Ontario 


The U. S. Lake Survey, in coordination with the Water Resources Branch, 
Department of Northern Affairs and National Resources, Canada, has deter- 
mined Lake Ontario end-of-month levels. From these data, the average 
monthly changes in Lake Ontario levels for the period 1934 to 1953 were 
computed and are tabulated in Column 2, Table 1. 

In the period of record, 1934-1953, the all-time mean annual low (1860- 
1957) of 243.54 feet was established in 1935, and the all-time annual high 
(1860-1957) of 248.01 feet was established in 1952. During this maximum 
secular climb of Lake Ontario, the average lake levels rose during the period 
January through May and fell during the period July through December; the 


TABLE 1 


Evaporation in Feet Computed by Water Budget Method 
Period of Record, 1934-1953 


Ey-2 = Sy-2 + R-2 + Ti-2 - %-2 + Pover water, _ 
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maximum monthly rises occurred in March and April, and the maximum 
monthly declines occurred in August and September. 


Runoff 


The runoff of precipitation from the drainage basin of Lake Ontario is an 
extremely important factor in analyzing lake level fluctuation. Preliminary 
studies have shown that there is a small lag from the time of precipitation 
until the runoff reaches the lake. . Although several of the tributary basins 
have a great number of small lakes and the basin of the Trent River is par- 
tially regulated, there does not appear to be any marked degree of storage in 
the Lake Ontario basin. Of the 34,795 square miles in the basin, the lake sur- 
face consists of 7,520 square miles and the land drainage area consists of 
27,275 square miles; therefore, it is not difficult to visualize why runoff is so 
important. 

The majority of runoff from the land area into Lake Ontario is in the form 
of stream flow from the major rivers tributary to the lake. The files of 
United States and Canadian agencies were examined for available stream flow 
records. Usable records for nine major rivers which include 16,500 square 
miles or 64 per cent of the Lake Ontario land area (excluding that area tribu- 
tary to the Niagara River above Niagara Falls) were found. Drainage areas 
with no available records were incorporated into the near drainage areas with 
surface records, see Figure 1. 


TRENT RIVER 


OSWEGATCHIE RIVER 


DUFFIN CREEK [7 
HUMBER RIVER 


CREDIT RIVER 


BLACK RIVER 


OSWEGO RIVER 


64%  GAGED AREA 
GENESEE RIVER 36% UNGAGED AREA 
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Figure 1. Drainage Areas, Lake Ontario 


The areas for which gaged records were not available were divided almost _ 
evenly between upland and coastal regions. Although there is not much 
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variation in the topography of the Lake Ontario Basin, except in the eastern 
portion where mountain foothills are located, care was taken to incorporate 
ungaged areas with the gaged areas that had the same general topographical 
characteristics. For instance, the ungaged area in the coastal region of 
northern Lake Ontario was incorporated with Duffin Creek which also lies 
within the coastal region. The ungaged upland regions east of the Moira 
River were incorporated with the Moira River basin. The size of the gaged 
| drainage areas of the nine rivers and the areas incorporated with the records 
of these rivers are given in Table 2. The Credit and Humber Rivers and Duf- 
fin Creek had limited records and the statistics of flow were mathematically 
extended to the period 1934-1953. 

The average monthly runoff for the lake basin measured in feet on Lake 
Ontario for the period 1934-1952 is given in Column 3, Table 3. 

Although there are certain to be some errors in measuring the total run- 
off, the gaged runoff is the most accurately measured variable in the hydrolog- 
ic cycle for, instead of being a point measurement within a large area (such 
as the measurement of precipitation), the gaged runoff effectively integrates 

_ the entire area from which the measured flow originates. 

Runoff is a complex phenomena which depends upon many variables. The 
observed monthly runoff for the Lake Ontario Basin is a function of the supply 
of water available for runoff and the losses which prevent the supply from 
reaching the lake. The supply of water available during a month is directly 
related to the precipitation during that month, the interflow, ground water, 
and channel storage, and the water equivalent of the snowpack; if any. The 
losses are due to interception, depression storage, evapotranspiration, and 
replenishment of soil moisture. 

Indexes of these most important variables affecting runoff were deter- | 
mined and a mathematical correlation of historical record with these indexes 
were obtained which gave results of rather high correlation. (1) 

The relation of runoff to precipitation proved to be a useful tool in gaining 
an insight into the runoff from the basin. This relationship for each of the 
drainage areas and for the entire basin is given in Table 3. Some very inter- 
esting results are shown. One, is the extremely large runoff factors for the 
months of March and April. This demonstrates conclusively that a knowledge 
of the storage of precipitation in the form of snow and the subsequent snow 
melt is of vital importance if a complete understanding of spring runoff is to 
be obtained. Another, is the magnitude of the yearly runoff factor for the 
Lake Ontario basin. The value of 48.0 per cent is appreciably higher than the 
value of 30.7 per cent obtained by John R. Freeman in his study, “Regulation 

of Great Lakes.” 


Precipitation over Lake Ontario 


The amount of precipitation over the lake affects the lake level immediately. 

Heretofore, the weighted precipitation catches recorded at selected weather 
stations on the perimeter of the lake have been regarded as the precipitation 
falling over the lake. The perimeter precipitation shows a marked reduction 
from the precipitation falling over the land drainage area, yet it seems that, 
due to the cooling effect of the lake in the summer and the fact that there are 
‘no diurnal heating or orographic effects, the precipitation falling over the lake 
is much less than indicated by the perimeter stations. 
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In late years there have been several concerted attempts to determine if 
the depth of precipitation falling over large bodies of water is of the same 
order of magnitude as the depth falling over the adjacent or surrounding land 
areas. In the Great Lakes area, precipitation has been recorded at St. James, 
Beaver Island, in northeastern Lake Michigan, for twenty-five years during 
the period 1911-1956. A comparison of the precipitation falling at Beaver 
Island with the precipitation falling on the adjacent land areas indicates that 
not only is the depth of precipitation less over Beaver Island, but that there 
is a seasonal variation. 


Northeastern Lake Michigan studies 


To investigate the theory that the depth of precipitation falling over large 
bodies of water is less than the depth falling over adjacent or surrounding 
land areas, a study of precipitation over northeastern Lake Michigan was 
instituted by the United States Lake Survey. The data for the study were ob- 
tained jointly with the U. S. Weather Bureau and with the cooperation of the 
U. S. Coast Guard. Storage precipitation gages were installed on six islands 
during November 1952, by the U. S. Weather Bureau, see Figure 2. These 
special island storage gages were inspected and the catch was measured 
semiannually. The results of the catch on the special island gages have been 
compared with a network of ten shore control stations. 

The topography of the six islands varies tremendously and the ground 
cover and exposure at the gage sites also differ greatly; yet, the catch at 


LAKE SUPERIOR 


ST. JAME! 
GULL I@ 9D 
BEAVER /SL, 


r\:) 
ILE AUX 
wud \ 


we 


$8. FOX iQ 


N. MANITOU ISL. 
8. MANITOU ISL.@ 


Legend: 


@ = PRECIPITATION 
STORAGE GAGES 


Figure 2, Precipitation Gage Locations, Northeastern Lake Michigan 
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ive of the six island stations for the period November 1952 to November 
11957, compares very favorably with that of St. James, Beaver Island. The 
erage precipitation of the five large island stations for the five year period 
is 103 per cent of the precipitation recorded at St. James. However, the 
catch at Ile aux Galets, the smallest of the six special island stations with 

an area of 2.9 acres, was appreciably less than the catch at the other island 
jstations. Ile aux Galets is a very small flat sand spit and is completely de- 
tate of vegetation, whereas the other island stations are fairly heavily vege- 
tated and relatively large in area. 

The greatest difference in precipitation between Ile aux Galets and the 
other island stations is in the winter period. This is believed to be caused 
‘by the extreme exposure of the gage, which reduces the amount of catch of 
snow. 
| The records during the summer period show that the depth of precipitation 
at Ile aux Galets is approximately the same as for South Fox Island. The 
average precipitation of Ile aux Galets and South Fox Island during the period 
May through October, 1952-1957, when comparable data are available, is 93 
per cent of that recorded at St. James. This is reasonable, because Beaver 
Island is large enough for definite diurnal heating to occur and the convective 
currents generated could cause additional precipitation as compared with the 
smaller islands. 

Nevertheless, to determine if the size and vegetative cover of the island 
had a definite effect on the amount of precipitation, an additional gage was 
installed in October, 1956, on Shoe Island, another sand spit 2.3 acres in size. 
Figure 3 shows the gage location on Shoe Island. The results of the first year 
of operation corroborate the reduced catch on Ile aux Galets, indicating that 
the rainfall over the small islands, and consequently over the lake, is re- 
duced from that recorded on Beaver Island. 

It is fully realized that five years and less are very brief periods on which 
to base specific conclusions regarding precipitation. But, general trends can 
be ascertained and approximate relationships can be determined. On the 
other hand, the twenty-five years of record at St. James are much more reli- 
able. The ratio of precipitation recorded at St. James to that recorded over 
the surrounding land area is given in Column 2, Table 4. These ratios show 
a definite reduction in rainfall, with the greatest reduction occurring in the 
Summer months. This is what may be expected, because there are fewer 
frontal storms in the summer months and those that pass over the lake are 
dampened due to the cooling effect of the lake surfaces. The convective type 
precipitation is pronounced over land during the summer period. In the win- 
ter, there are more frequent passages of frontal storms, no dampening effect 
by the lake, and no squall-type storms. Therefore, one would expect the 
precipitation over the lake to approximate that occurring over the land. The 
difference in the air and water temperature in the vicinity of Beaver Island 
is illustrated in Figure 4. It can be seen that the average monthly air tem- 
perature is higher than the average monthly water temperature during the 
period April through August, accounting in part for the reduced rainfall over 
the lake during this period. 

In a previous paragraph it was stated that the precipitation over the lake 
was indicated to be less than that recorded at St. James. From the records 
of South Fox Island, Ile aux Galets, and Shoe Island, it seemed reasonable and 
conservative to assume the reductions shown in Column 3, Table 4. Thus, it 
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Figure 3. Precipitation Gage, Shoe Island 


was assumed that the product of the two ratios, the ratio of St. James to land 
area and the ratio of small island stations to St. James, would give an approx 
mate solution to the reduction in precipitation falling over northeastern Lake 
Michigan. A plot of a smoothed curve drawn through the values of the produc 
of the two ratios, or the ratio of precipitation over water to precipitation ove. 
land, appears on Figure 5, and the “smoothed” monthly values used to reduce 
land precipitation over northeastern Lake Michigan are given in Column 5, 
Table 4. 

In his paper entitled, “The Influence of Lakes and Urban Areas on Radar 
Observed Precipitation Echoes,” Bulletin of American Meteorological 
Society, J. E. Pearson gives the conclusion that, “Observation of radar film 
of echoes in the vicinity of or over Lake Michigan revealed evidence that the 
lake, in general, discourages the formation of air mass showers.” This fact 
is definitely evidenced by the outline of radar echoes as depicted in Figure 6. 
In addition, data collected by the Jacksonville District, Corps of Engineers, < 
precipitation stations on Lake Okeechobee, the second largest fresh waer la 
in the United States with an area of 730 square miles, reveal a reduction in 
rainfall over the lake compared with the catches on the perimeter of the lake 
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TABLE 4 


Precipitation Ratios, Northeastern Lake Michigan 


2 4 5 

St. James to | g Small Islands | "Smoothed" Small 

Surrounding Rae es to Surrounding Islands to 

Land ; Land Surrounding Land 
January 25 -90 -90 89 
February . 86 86 285 
March Aids «76 oth) 
| Borin .80 76 .76 
May 81 5th sie) 
| June -85 279 ~74 
July 3 hee .69 .69 
August ays .60 62 
September .80 aie a2, 
October ~ 96 .89 85 


ofS -90 


November oe 


December .88 88 90 


Yearly 
Average 


*Hypothetical 


[The U. S. Weather Bureau office in Key West has approximated their rainfall 
is 60 per cent of that falling on the Florida mainland. Thus, there is conclu- 
sive evidence that the precipitation falling over large bodies of water is re- 
luced. 4 

_ One may reasonably ask if the ratios obtained for northeastern Lake 
Michigan are applicable to Lake Ontario. Evidence indicated that approxi- 
nately the same ratios are applicable. Generally, the same prevailing winds 
iold for the whole Great Lakes region. A comparison of precipitation at the 
erimeter stations with the land drainage area stations of Lake Michigan and 
| 


} 
| 
24 February, 1959 HY 2. 


Tw-Ta, °F 


MONTHS 


Legend: 


Tw- AVG. LAKE ONTARIO 
rata TaA- MEAN TORONTO,OSWEGO,TRENTON 1937-56 


© Tw- AVG. LAKE MICHIGAN, VICINITY OF ST. JAMES ISL 
~ \ TA- MEAN ST. JAMES, PERIOD OF RECORD 


Figure 4. Average Monthly Air-Water Temperature Differentials 
in Degrees Fahrenheit, Lakes Michigan and Ontario 


Lake Ontario showed that the ratios are remarkably the same. As shown in 
Figure 7, the reduction in precipitation over both lakes follows the same pat- 
tern. Furthermore, a comparison of the average monthly water temperature 
minus the average monthly air temperature revealed that almost exactly the 
same pattern is followed by both Lakes Michigan and Ontario, see Figure 4, 
above. This indicated that the lakes discourage the formation of air mass 
showers during the same period of the year. Accordingly, factors developed 
for reducing the amount of precipitation falling over northeastern Lake 
Michigan will apply to Lake Ontario. 

The average monthly precipitation falling over the Lake Ontario land 
drainage area is constant, varying but slightly from the average monthly 
mean of 0.24 foot. The average monthly over-land precipitation is shown in 
Column 5, Table 1, above. This average monthly precipitation over land, 
multiplied by the ratios previously computed for reducing the precipitation 
over land to the precipitation over water Column 5, Table 4, gave the aver- 
age monthly depth of precipitation falling on Lake Ontario, see Column 6, 
Table 1. 


Inflow to Lake Ontario from Lake Erie 


The inflow to Lake Ontario from Lake Erie is the flow of the Niagara 
River plus the diversion of water through the Welland Canal in Canada. The 


2 LAKE ONTARIO 25 


1S THAT OF NE LAKE 
ORAINAGE AREA. 


WATER PRECIPITATION IS THAT RECORDED 
AT ST. JAMES, BEAVER ISLAND., 25 YEARS I9II-56, 
ADJUSTED FOR RAINFALL OVER SMALL ISLANDS. 


lod Fo Ma) Ate ESOS AC Sate Oo ND 
MONTHS 


Figure 5. Ratios, Precipitation Over Water to Precipitation Over ‘Land 


Welland Canal diverts water for navigation and power from Lake Erie at Port 
Colborne and returns it to Lake Ontario at Port Weller and Port Dalhousie. 
The Niagara River flows used in this study were computed by the new Morrison 
Street rating curve. The Morrison Street gage is located in the Maid of the 
Mist Pool below Niagara Falls and was the only gage on the Niagara River 
which was comparatively free from ice and weed effects during the period of 
record. It is interesting that the flows used in this study differ from the pub- 
lished flows of the Niagara River which occurred at high Lake Erie levels by 
4s much as 10,000 cfs per month. Revision of the published Niagara River 
“lows is being studied by the United States Lake Survey. 


Sutflow of the St. Lawrence River 


_ The Coordinating Committee on Great Lakes Basic Hydraulic and 
Tydrologic Data(2) has accepted the premise that the Rapide Plat weir is the 
est available meter of the flow of the St. Lawrence River and that the flow 
ver the weir can be correctly measured by utilizing the stage at the Lock 25 
rage. The committee, therefore, adopted the Lock 25 stage-discharge rela- 
ionship as the basic relationship for the determination of the river flow. The 
sake Ontario outflows utilized in this study are those adopted by the Coordi- 
ating Committee on Great Lakes Basic Hydraulic and Hydrologic Data. 
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Figure 6. Outlines of Precipitation Echoes Taken at 30-Minute 
Intervals. Courtesy of Illinois State Water Survey 
Division, Urbana, Illinois 


The average monthly inflows to Lake Ontario from Lake Erie minus the 
average monthly outflows of the St. Lawrence River, in feet of Lake Ontario 
levels, for the period 1934-1953, are shown in Column 4, Table 1, above. 
These values do not vary greatly and there is excellent correlation between 
the inflows and outflows when they are analyzed statistically. 


Other factors 


All other factors which affect the amount of evaporation computed by a 
water budget equation not previously discussed were included in the item 


U,-2. The most important of these factors were the ground water flow and 
the thermal expansion of water. 
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Figure 7. Ratio, Precipitation at Lake Perimeter Stations to 
Precipitation Over Land, Lakes Michigan and Ontario 


sround water flow 


The average level of Lake Erie is 572.3 feet, whereas, the average level of 
sake Ontario is 246.0 feet. It was possible that there might be a considerable 
inderground flow between the two lakes due to the very large potential head. 
This matter was discussed with members of the U. S. Geological Survey and 
he consensus was that there was not an appreciable ground water flow into 
sake Ontario. 


*hermal expansion of water 


As a mass of water is heated and cooled, it will expand and contract. The 
reater the depth, the larger the thermal expansion and contraction. The 
verage depth of Lake Ontario is 264 feet. However, as already shown, a 
efinite thermocline exists in Lake Ontario. This thermocline varies in depth 
oth with time and position, but the maximum depth does not exceed 50 feet, 
slow which the water is generally at a temperature of 399 F. The maximum 
lonthly average water temperature difference between two consecutive 
ionths was 11° F, see Table 4 above. The lake expands in a vertical direc- 
on only, and the expansion is a function of the change in average density of 
ie water above the thermocline. The maximum amount of expansion and con- 
action is about .01 foot per month. Consequently, the thermal expansion of 
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water was not considered in the water budget equation. 

In addition to the ground water flow and the thermal expansion of water, 
consideration was given to other probable factors which might not be negligi 
ble. These factors were grouped under Uj_9 in Equation (1) and were con- 
sidered to be zero. 


Water budget evaporation 


The average monthly evaporation of Lake Ontario for the period 1934 to 
1953 computed by the water budget method was obtained by substituting the 
values of the normal change in levels, runoff, inflow minus outflow, and pre- 
cipitation over the lake in Equation (1). The values so obtained are given in 
Column 7, Table 1 above. Note that the minimum average monthly evapora- 
tion occurred in June and the maximum average monthly evaporation oc- 
curred in October. The reason for this will be explained below. The total 
annual average evaporation for the period was 2.21 feet. 


Mass Transfer Method 


The velocity of a viscuous fluid at the zone of contact with a solid bounda: 
must have the same velocity as the solid boundary. Thus, a free stream of 
air moving with the velocity U past a fixed solid boundary must have a veloc 
' ty of zero at the boundary and a velocity of U some distance away. Where a 
velocity gradient is found, the turbulent motion leads to a transport of mome 
tum across surfaces normal to the velocity gradient. A shearing stress is 
exerted on a surface across which a transfer of momentum takes place. Th 
relatively thin layer where the velocity gradient is large is called the bound 
ary layer. Not only is the character of the fluid flow, laminar or turbulent, 
important, but the nature of the boundary has an important function. The flo 
of wind over Lake Ontario was considered to be turbulent and the surface 
aerodynamically rough. 

The general equation for shearing stress is: 


bet 
read (2 


where 
Y is the shearing stress. 
A is the coefficient of turbulent interchange of momentum. 


a is the gradient of horizontal wind velocity. 


The transfer of water vapor through the atmosphere was considered to o1 
cur in the same manner as the transfer of momentum. Thus: 


Tie eh get (3 
where 
E is the evaporation. 


A is the coefficient of turbulent interchange of moisture. 
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ays : 
io is the gradient of moisture concentration. 


The specific humidity is defined as 


‘4 622 
QS eee (4) 


Pa = .378e 
Mibere 


e is the vapor pressure of the air 
Pa is the atmospheric pressure. 


ince Pa is of the order of one hundred times e, there is little loss of ac- 
uracy by substituting 


_ .622e 
Pa (5) 


Since the coefficient of turbulent interchange of moisture was considered 
to be equal to the coefficient of turbulent interchange of momentum, the 
evaporation was considered to be a function of the friction velocity, Ux, and 
the difference between the vapor pressure of the air and the vapor pressure 
of saturated air at the water surface. 


He= Baby Ave (6) 


where Ux is the friction velocity and is equal to {7 / 2 


Equations for the value of the shear stress of the wind over water under 
all conditions of temperature stability have been determined.(3) In using 
these equations, it was assumed that when the water surface is covered with 
waves different equations must be used than when there are no waves on the 
surface. There appeared to be a critical value of U, for which the water 
surface becomes unstable, and waves form. 

Excellent data obtained from water-loss investigations at Lake Hefner 
were used to correlate measured evaporation with that computed by Equation 
(6) above. The equation developed from the Lake Hefner data for evaporation 
in inches per day which was considered to be applicable to Lake Ontario was 


BE =7,000733 U,.(e,-- 6. ) (7) 


where 
U, is in centimeters/second. 


e is in millibars and the vapor pressure of the air is measured at 2 
meters elevation. 


The correlation factor of Lake Hefner measured evaporation to that com- 
uted by the above formula was 0.90. The equation also appeared to give good 
‘esults for the Lake Mead studies. 

- It is an extremely difficult task to obtain wind and temperature data for 
Lake Ontario. The friction velocity (Ux) is a function of the wind velocity over 
yater at a given height and the atmospheric stability. The measure of atmos- 


heric stability used was the difference between the air and water surface 
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temperatures (TA - Tw). There are only two anemometers operating the 
year round near Lake Ontario where the wind velocities measured over land 
can be adequately reduced to over-lake wind velocities. Both of these ane- 
mometer locations, Toronto and Cobourg, are in Canada, and, unfortunately, 
the Cobourg anemometer was discontinued in 1952. The recorded wind 
speeds were first reduced to eight meter level wind speeds and then were in 
creased to over-water winds in accordance with experience gained in relatir 
vessel weather observations of Great Lakes anemometer -equipped vessels 
with land stations.(4) The factors by which the over-land wind speeds were 
increased depended upon the atmospheric stability; these are given in Table 
5. 


TABLE 5 


U 
Ratio of Over-Water to Over-Land Wind Speeds, jaa 
L 


The vapor pressure differential (Ae) is a function of the temperatures of 
the air, the relative humidity of the air, and the water surface temperature. 
The air temperatures and the dew points were derived by arithmetically 
averaging the values recorded at Toronto, Trenton, and Oswego, for the 
period 1937 through 1952. The water surface temperatures were derived 
from the work of Freeman, Millar, and the U. S. Lake Survey, see Table 6; 

The evaporation computed by the mass transfer equation, Equation 7, is 
given in Table 7. 

Using the mass transfer method, the minimum average monthly evapora 
tion occurred in May and June and the maximum average monthly evaporati 
occurred in September. The difference in vapor pressures in May and June 
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TABLE 6 


Lake Ontario Water Surface Temperatures 
Degrees Fahrenheit 


U. S. Lake Survey? Adjusted 
Temperatures 


aN 


Freeman Millar@ 


January 

February 32 35 
March 33 35 
April 39 BY. 
May 46 42 
June Die SPA 
July 67 


August 


September 


October 


November 


December 


dipegulation of the Great Lakes," report by John R. Freeman. 


2uSurface Temperatures of the Great Lakes," F. Graham Millar, Journal 
of the Fisheries Research Board of Canada, December 1952. 


3 Archives, U. S. Lake Survey, Detroit, Michigan 


shows that there was no appreciable gradient of moisture concentration and, 
consequently, no net loss of water mass. In fact, in some years, there was a 
negative evaporation in the months of May and June. 


Discussion 


; 

: All data from known sources were used in this study. However, the avail- 
ole data were not of the highest order, consequently, the computed evapora- 
on is subject to error. New installations must be established to collect 
nore complete and accurate meteorological and hydrological data in the 
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TABLE 7 


Evaporation in Feet Computed by Mass Transfer Equation 


Over Water 


Boas 8 anspor tia Sli ain eae saaeage roct/icil 
Centimeters/second 
January 832 43.4 081 
February 806 42.7 -077 
March 636 34.1 -073 
April 588 36.0 
May 398 26.9 
June 338 22.4 .013 
July 374 21.7 
August 392 21.6 olor 
September 455 22.8 
October 479 23.8 .150 
November 605 3h. 2 - 100 
December 818 2 7 ie 
Yearly 
Total 2.00 


Great Lakes area if the most beneficial regulation plans are to be effected. 

A comparison of the evaporation computed by the water budget and the 
mass transfer methods is given in Table 8. 

The results compare rather favorably and the trend of evaporation is the 
same, although there are minor differences in the amount of net evaporation. 
The evaporation in the spring months has been shown to be negligible becaus 
there is no appreciable gradient of moisture concentration. The maximum 
evaporation occurs in the fall. 

There is definite evidence that the precipitation falling over the surface o 
the Great Lakes is considerably less than the precipitation falling on the ad- 
jacent areas. The fact that the precipitation over Lake Ontario is some 20 
per cent less than the precipitation recorded at land stations throughout the 
basin means that the values of evaporation are reduced in the same amount 
when computed by the water budget equation. As a result, the annual amount 
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TABLE & 


A Comparison of Computed Evaporation, in Feet 


Water Budget Mass Transfer 
Method Method 


[ees 
January 19 523} 


February -10 Rae, 
March -06 6) 
April 04 .09 
May : 04 .00 
June -00 .02 
July a .18 


August .28 +22 


September ere rete 


October 38 AE: 
November ne) AP) 


December wee 


of net evaporation from Lake Ontario computed by both the water budget and 
the mass transfer methods is of the order of 2.0 feet. 


REFERENCES 


1. Major Ira A. Hunt, Jr., A Technique for Forecasting Observed Runoff from 

Large Areas, U. S. Lake Survey, (Detroit, 1958). 

2. Appointed by the U. S. and Can. governments to develop Great Lakes basic 

_ hydraulic and hydrologic data and to establish the basis for future develop- 

ment of such data by these agencies. 

Capt. Ira A. Hunt, Jr., Effets du vent sur les nappes liquides, (The Effects 

of the Wind on Liquid Surfaces), Univ. of Grenoble, (Grenoble, France, 

+1954). 

4, Major Ira A. Hunt, Jr., Relationships Between Over-Water and Over-Land 
Wind: Lake Erie, U. S. Lake Survey, (Detroit, 1958). 


oe 


‘ at a YS 


* 
o> 


> 


1935 February, 1959 HY 2 


EE 


Journal of the 
HYDRAULICS DIVISION 
Proceedings of the American Society of Civil Engineers 
—— 


PROBLEMS CONCERNING USE OF LOW HEAD RADIAL GATES 


Thomas J. Rhone,! A.M. ASCE 


ABSTRACT 


Three of the principal hydraulic features of low head radial gates are dis- 
cussed. These are (1) a general discharge equation, (2) the effect of gate seat 
location on discharge capacity and pressure distribution along the spillway 
surface, and (3) side and bottom seals. 


INTRODUCTION 


The radial-type gate originated in France about 100 years ago.(1) The 
earliest recorded use was by the French Engineer Poiree on the Seine River 
in 1853. In the 1860’s, another French Engineer, Mongel Bey, used a cast- 
iron, radial-type gate in the Delta Barrage on the Rosetta Branch of the Nile 
River. The Poiree gates were 28.7 feet wide by 3.3 feet high; the Delta 
Barrage gates were 16.4 feet wide by 16.7 feet high with the concave face 
turned toward the reservoir. 

The radial-type gate was patented in the United States in 1886. Since then, 
this type of gate has become widely used and has increased in size. Three 
114 feet wide by 26 feet high gates are used on Horseshoe Dam Spillway near 
Phoenix, Arizona, and it is proposed to use four 40-foot wide by 52.5-foot high 
radial gates for reservoir control at the new Glen Canyon Dam. 

Radial or Tainter gates are probably the most widely used crest control 
gates. They are particularly well adapted to crest control because of their 
Simplicity of design, construction, and installation. The overall economy and 
efficiency of this type of gate result from the radial bearing. The thrust from 
the waterload is carried to two trunnion bearings where it offers little resist- 
ance to the operation of the gate. This action can be compared to the hub of a 
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wheel where all of the force on the rim is transmitted to the center. In lift- 
ing the gate, there is some force to be overcome at the trunnion, and the side 
seals may offer considerable resistance in the form of friction. However, 
the opening of a radial gate requires less hoist capacity than a slide gate with 
face bearing. In addition, radial gates are more adaptable for automatic con- 
trol apparatus. 

Another desirable feature of a radial gate is that it needs no gate slots. 
High head flow past indented slots has in some installations produced cavita- 
tion damage to the pier and spillway surfaces. One instance where gate slots 
were the apparent cause of cavitation to spillway and pier side walls is 
shown in Figure 1; several other cases have been reported in technical 
literature.(2) Preventive action has included the use of gate followers to fill 
the slots at small gate openings, and the use of offsets in the pier side walls 
downstream from the slots. 

Despite the apparent simplicity and wide use of radial gates, very little 
seems to be known about their operating characteristics. Specifically, there 
has been very. little information published on the discharge capacity, the ef- 
fect of gate location on pressure distribution on the spillway face, downpull 
produced by flow under the gate, and the type of seals to be used, and other 
hydraulic problems. 

The purpose of this paper is to discuss some of the hydraulic features of 
radial gates, including the side and bottom seals; to present hydraulic data on 
some representative gate installations; and to create interest in the remain- 
ing problems, with the hope that additional material will be presented by the 
readers of this paper. Further, it is hoped that the tests described will 
stimulate hydraulic research on radial gates by those in a position to do 
hydraulic model or prototype investigations. 

Material in this paper has been limited to the so-called low head radial 
gates used to control flow over spillways. Top seal radial gates, where the 
normal operating head greatly exceeds the gate height, and the small gates 
used to control canal flow, are not covered. 


Discharge Characteristics 


Discharge Determination 


Although the radial gate has been used for regulating flow for over a 
hundred years, a general method for predicting exact flow quantities has not 
been derived except for one type of installation. When the gate is used in a 
flat-bottom, rectangular section and the discharging jet is supported by the 
flat floor, the flow characteristics are known both from mathematical and 
experimental analyses. The discharge coefficients derived from these studie: 
have been confirmed by several experimenters and are accurate for their 
specific purpose.(3,4) 

Many specific discharge determinations are on record; usually, when a 
spillway has been model-tested, calibration curves for both free and gate- 
controlled flow have been obtained. From an analysis of these ‘calibrations, ~ 
the free flow discharge coefficients for almost any shape of overflow section 
can be determined.(5) However, a general equation which will provide the 
exact discharges for radial-gate-controlled flow is still unpublished. 

The lack of a general equation can probably be explained by the many vari-. 
ables that affect the flow pattern. Some of these are in the spillway approach 
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and include channel width and depth. Other variables include the shape of the 
spillway profile; the geometry of the gate, including the radius, width, height, 
trunnion location, and the location of the gate seat; the location, spacing, and 
type of piers; and the method of determining gate opening and reservoir ele- 
vation. The combinations of these variables are almost countless, and the 
effect of any one of them is difficult to determine from a limited number of 
investigations. 

A brief review of four methods which have been used to determine the flow 
quantities under radial gates indicates the different approaches used in at- 
tempts to obtain a general solution. It should be pointed out, however, that 
each provides only approximate discharges; if an accurate determination is 
necessary, either model studies or field calibrations should be made. 


: aay ey fea 3/2 e: 
The first equation is Q = 2/3 CL V2g h, <-hy where “C” is the 


coefficient of discharge listed as a function of gate opening and reservoir 
head. The definition of the symbols and suggested values for “C” are shown 
in Figure 2. From the graph, it is apparent that for h;/d values less than 
about 2.2 there is a wide range where “C” can vary by as much as 20 percent. 
The coefficient curve shown as the heavy line was derived from hydraulic 
model studies of the spillways listed on Figure 2. The different crest shapes 
and gate arrangements probably account for some of the spread. However, it 
is known that where the head on an orifice is small in comparison to the 
height of the orifice, there is an appreciable difference between the discharge 
obtained by using an average head and the discharge obtained by taking into 
consideration the variation in head. Depending on the actual values of the 
orifice height and head, the differences in discharge may be large or small. 
Thus, the region below hj/d = 2.2 is a transition region, a region notoriously 
difficult to evaluate in all hydraulic flow problems. 

A second method of determining discharge quantities takes into considera- 
tion the angle 6 of lip of gate to horizontal and the angle A described by the 
intersection of a horizontal line through the pivotal point of the gate and the 
radius of the gate drawn from the same point to the center of the gate when it 
is in a closed position, Figures 3 and 4. The value of the coefficient of dis- 
charge for controlled flow is related to the known coefficient of discharge for 
free flow for different values of 9 and ratios of gate openings d to depth of 
flow hy. Plots for two different values of @ are shown in Figures 3 and 4. 
The rate of flow is then obtained by applying the formula 


9/3 eeniahaey 
Q = CL/hy —ho 


where hy = vertical distance from crest to water surface 
hg = vertical distance from bottom of gate to water surface 
L 


The coefficients were derived from the results of many model studies. 
The third method is similar to the second in that it relates the controlled - 
flow discharge coefficient to the free flow coefficient, Figure 5. The effect ot 
gate radius, trunnion location, and gate opening is also recognized as having 
an effect on the coefficient and is evaluated by the parameter _1 . The 
sin 6 
values for “C” were obtained from analyses of calibrations of model spillway 
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DISCHARGE UNDER A RADIAL GATE 
FIGURE 5 


The discharge formula for this method is: 


1 oi omens 
Q = CL sine Ay -hy Hi 


Values for CL and 1/sin @ are given in Figure 5. : 

The fourth method, Figure 6, takes into consideration the crest shape, lo- 
cation of gate seat, the angle 6 formed by the tangent to the gate lip and the 
tangent to the crest curve at the nearest point of the crest curve, and the net 
gate opening which is the shortest distance from the gate lip to the crest _ 
curve. From the geometry of the gate the angle 6 is determined and from tt 
curves in Figure 6, a coefficient of discharge is determined. These values 


8 IN DEGREES 


RADIAL GATES 


*0.60 0.65 0.70 0.75 
DISCHARGE COEFFICIENT-C 


EQUATION FOR DISCHARGE 


Q=sC DEV Zone 


D = Net gate opening 

L = Crest width 

H = Head to center of gate opening 

For C use dashed line when gate 
seats on crest and solid line when 
gate seats below crest. 


‘DISCHARGE UNDER A RADIAL GATE 


FIGURE 6 


43 


: 
44 February, 1959 HY 2 


for “C” were also obtained from model investigations. The discharge is then 
computed from the formula: 


Q = cDL ¥2eH 


where Q = discharge in cfs 

C = coefficient of discharge 

D = net gate opening 

L = crest width 

H = head to the center of the opening 


Table 1 shows Canyon Ferry Dam Spillway discharges computed from the 
four methods described, compared with the discharges obtained from model 
and prototype calibration. The Canyon Ferry Dam Spillway is a radial-gate- 
controlled spillway located on the Missouri River about 17 miles east of 
Helena, Montana. The spillway was calibrated by model studies in 1953, and 
field discharge measurements were obtained in 1956.(6) 

The first and fourth methods give 5 to 11 percent higher discharges than 
actually measured. The second and third methods give discharges from 
2 percent over to 7.5 percent under the measured values. The small gate 
opening shows a greater variation than the larger openings. 

The obvious conclusion that can be drawn from Table 1 is that a satisfac- 
tory method for computing the discharge through partially open radial gates 
has not yet been derived. A method of computing discharges is required whic 
will be in better agreement with model tests and field calibrations. Examina- 
tion of the four discharge equations shows that in three of them the discharge 
varies as hy3 2 and in the fourth as H!1/2, A more nearly correct method 
might be to vary this exponent with changes in h,/d. The change from free 
weir flow to orifice flow, or vice versa, would, thereby, be taken into ac- 
count. The transition range occurs when hj/d is less than 2.2, Figure 2. 
When the ratio is greater than 2.2, true orifice flow occurs. Thus, in the 
free flow range the exponent for the head in the discharge formula is three- 
halves; in the transition range, the exponent should vary between three-halve 
and one-half; for orifice flow, the exponent is one-half. 


Effect of Gate Position on Discharge Coefficient 


In 1953, hydraulic model investigations were made to determine the effect 
on pressure distribution along the spillway profile caused by changing the lo- 
cation of the gate seat.(7) 

For this study, 5 crests were investigated. With each crest, 4 or 5 gate 
seat locations were tested. The locations were (1) on the crest, (2) and (3) 6 
inches vertically below the crest both upstream and downstream from the 


3 ft : 32.32: h,hoo: 4,450 : 4,870 : 4,110 : 4,060 : 4,880 
6 ft : 32.00 : 8,450 : 8,400 : 9,240 : 8,370: 8,270: 8,9h0 
10 ft : 32.45 : 13,800 ; 13,750 : 14,720 : 14,070 : 13,980 : 14,510 


*Gate openings and heads are measured vertically above gate 
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crest axis, (4) 12 inches below the crest on the downstream side, and (5) the 
design gate seat location when it did not correspond to one of the other four. 

The five crests and the model scale ratio at which each was studied are 
(1) Ross Dam 1:24.14, (2) Altus Dam 1:20.70, (3) American Falls Dam 1:18.35, 
(4) Bhakra Dam 1:50.56, and (5) Canyon Ferry Dam 1:42.92. The crest pro- 
files are shown on Figures 8 to 12. 

As a part of these tests, data necessary to compute the coefficient of dis- 
charge for several gate openings with each gate location were obtained. The 
coefficients were computed using the first method described. It was found 
that when the gate seats either upstream from the crest or on the crest, the 
coefficient of discharge is higher than when the gate seats downstream from 
the crest. This was true for all crests; in some cases, the difference was 
as much as 15 percent. 

Some of the difference might be explained by the fact that the head and gate 
Opening were measured from the gate seat instead of the crest. However, the 
head and gate opening values are the same whether the gate seats 6 inches be- 
low the crest in an upstream or downstream direction from the crest. In all 


_ cases, the coefficients for the upstream position were always higher. An ex- 
ample of these coefficients is shown in Figure 7. 


Conclusions and Recommendations 


Four different discharge equations have been used to compute the flow 
quantities under the radial gates of the Canyon Ferry Dam spillway. The 
discharge coefficients used in each equation were based on different model 
studies; each equation placed emphasis on different geometric features of the 
structure. In one equation, the head was raised to the one-half power; in the 
others, the head was raised to the three-halves power. The flow quantities 
obtained from these methods were compared with actual model and prototype 
measurements; the computed flows were found to vary from the measured 
quantities by as much as 11 percent. 

All of the methods described would probably give discharges sufficiently 
accurate for design purposes or for discharge determinations where water 
has little value and a rough estimate of the flow would suffice. However, when 
it is necessary to have an accurate discharge determination for flood routing 
or to measure valuable irrigation water, as in the western United States, 
either a model or prototype calibration is necessary. 

An accurate formula could be derived by combining theory and data analy- 
sis from coordinated hydraulic investigations on model or prototypes. The 
amount of existing data available from hydraulic laboratories might provide 
sufficient material for a preliminary analysis and if these data were supple- 
mented by models specifically designed for discharge studies, the effect of 
each variable could probably be evaluated. Such models would have to be 
sufficiently large that changes in the geometry of the gates, piers, and over- 
fall section would make measurable changes in the flow. The crest length 
should be sufficiently long to allow installation of several piers, and provision 
should be made to test several crest shapes, pier shapes, and gate sizes. It 


‘would be desirable to incorporate an adjustable weir in the model headbay, 
‘which would maintain a constant reservoir level to facilitate the flow measure- 
‘ment in the transition range between free flow and orifice flow. Standard 
‘measuring devices for determining gate openings, discharge quantities, and 
heads on the crest would be adequate. With such a model, the effect of any of 
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seat 2.43 ft. below crest,on downstream side 
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the variables such as approach depth, pier shape, or trunnion location could 
be evaluated. 


Pressure Characteristics 


A standard spillway crest is usually designed to conform to the shape of 
the lower nappe downstream from a sharp-crested weir. This profile is 
called the datum shape.(5) When gates are to be used on the spillway, the 
Same crest shape is used from the upstream face to the gate seat. Down- 
stream from the gate seat, the crest profile is sometimes made to corres- 
pond to the theoretical trajectory of a jet flowing from a 1-foot high rectangu- 
lar orifice at the design head. For the same design head, the crest designed 
for free flow will be steeper than a crest designed for a 1-foot gate opening. 

Several factors other than hydraulic phenomena are usually considered 
when determining the location of the gate seat. If the gate seats on the crest, 
the theoretical crest curve will be comparatively flat and the spillway will 
contain considerably more concrete than is usually needed for stability. If 

the gate seats below the crest, the profile will be steeper, but a higher gate 
is necessary to retain the same storage pool. Structural and architectural 
considerations sometimes dictate the gate seat location. If a bridge spans 
the spillway, it might be necessary to place the gate downstream from the 
crest to provide space for the gate lifting mechanism or to provide clearance 
when the gate is raised. The gate might have to be located downstream to in- 
sure that the trunnion will be above the water surface. 

Many hydraulic designers believe that when the gate is seated on the crest 
there is a tendency for the jet to spring free of the crest at high reservoir 
elevations and small gate openings; this may result in subatmospheric pres- 
sures and possible cavitation damage to the spillway. The designers also be- 
lieve that placing the gate seat downstream from the crest results in the flow 
being directed downward; the tendency of the jet to spring free of the crest is 
lessened. However, observations of prototype spillways indicate that cavita- 
tion damage to the spillway surface occurs downstream from the region where 
the control gates have any effect on the flow. No damage to spillway crests in 
the vicinity of the gates has been found which can be attributed to the crest 
shape and gate seat location. 

-A realistic solution to these problems and conflicting requirements is 
needed. The following discussion may help to clarify the problems. 


Effect of Gate Location on Pressure Distribution Along 
a Spillway Profile 


The laboratory investigations used to determine the effect of gate location 
on discharge coefficients were also used to determine the effect of the gate 
position on the pressure distribution along the spillway surface.(7) For the 
pressure tests, 10 piezometers were installed along the crest profile. One 
piezometer was placed at each tested gate seat and the others were equally 
spaced downstream along the spillway profile. The five crests used in the 
investigation provided a wide range of steepness and gate seat locations. 

For convenience, each of the crests will be discussed separately giving a 
description of the crest and the results of the tests. The crest profiles, 
datum shapes, piezometer locations, and pressure profiles for representa- 
tive gate openings are shown on Figures 8, 9, 10, 11, and 12. 
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Ross Dam, Figure 8 


The design gate seat is 4 inches below the crest. The datum shape is 
slightly steeper than the design shape. A comparison of the design shape and 
datum shape indicates that the overfall section from the gate seat downstream 
was designed from the theoretical trajectory of a jet from an orifice. From 
the reasoning previously outlined, pressures along the profile should increase 
when the gate seat is downstream from the design location and reduce when 
the gate seat is upstream. 

The pressure profiles on Figure 8 show that for the same gate openings 
and heads, the pressures for all gate seat locations downstream from the de- 
sign gate seat are practically the same. However, when the gate is seated 
upstream from the crest, the piezometers from the crest down to the second 
piezometer below the design gate seat showed pressures lower than were ob- 
tained for the other gate seat locations. Another trend which may be ob- 
served from the pressure profiles is that the zone of subatmospheric pres- 
sures moved upstream as the gate seat moved upstream. However, for all 
gate positions, the downstream piezometer (No. 9) showed pressures near 
atmospheric. 


Altus Dam, Figure 9 

The design gate seat is on the crest. The datum shape is the same as the 
design shape, indicating that the profile was obtained from. the lower nappe 
shape of a sheet of water passing over a sharp-crested weir. According to 
reasoning, the pressures along the profile should have the highest values 
when the gate seat is 12 inches below the crest and should become smaller 
as the seat is moved upstream. The pressure profiles in Figure 9 show this 
to be true in that for both downstream gate seat locations the piezometers 
downstream from the gate seat indicate nearly atmospheric pressures for all 
gate openings. Also, when the gate seats upstream from the crest, the piezo- 
meters from the gate seat to Piezometer No. 6 indicate subatmospheric pres 
sures. The piezometers downstream from No. 6 all show near atmospheric 
pressures. Generally speaking, Piezometers No. 6 through 9 showed nearly 
atmospheric pressures for all gate seat locations at all gate openings. Piezo 
meters No. 2 through 5 were affected by the gate seat location and the pres- 
sure was reduced at each piezometer as the gate seat was moved upstream. 


American Falls Dam, Figure 10 

The design gate seat is on the crest. The datum shape is considerably 
steeper than the design shape, indicating that the profile from the crest down 
stream was obtained from the theoretical trajectory of a jet flowing from an 
orifice. The pressures should be near atmospheric for small gate openings 
when the gate is at the design location, for downstream gate seat locations 
the pressures should be higher, and for upstream locations the pressure 
should be lower. The pressure profiles on Figure 10 indicate that perform- 
ance agrees with the reasoning. The only pressures significantly subatmos- 


pheric were those at the crest piezometer when the gate seat was upstream 
from the crest. 


Bhakra Dam, Figure 11 

The design gate seat is 2.43 feet below the crest on the downstream side. 
The datum shape is steeper than design shape. For this crest, the highest 
pressures should occur with the gate at the design location and they should . 
become successively lower as the seat location is moved upstream. Also, 
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ALTUS DAM CREST 
FIGURE 9 


the piezometers should indicate near atmospheric pressures when the gate 
seats on the crest. The pressure profiles on Figure 11 show this to be es- 
sentially true. The piezometers downstream from the gate seat show slight 
reductions in pressure as the seat is moved upstream and when the gate seats 


on the crest and upstream from the crest, the pressures are a little above an 
below atmospheric. 


Canyon Ferry Dam, Figure 12 3 

The design gate seat is 6 inches below the crest. The datum shape is 
steeper than the design shape. The relative shape of the two profiles indi- 
cates that the profile downstream from the design gate seat was derived from 
the theoretical trajectory of a jet flowing from an orifice. Therefore, the 
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(2) GATE SEAT AT PIEZOMETER 4, 12" BELOW CREST Gate open 14.3% of design head 
Gate open 428 % of design head 
Gate open 57.2% of design head 


EFFECT OF GATE POSITION ON PRESSURE DISTRIBUTION 
BHAKRA DAM CREST 
FIGURE II 


pressures downstream from the gate should be near atmospheric when the ~ 
crest is at the design location, higher when the gate seat is downstream and 
lower as the gate seat is Moved upstream. The pressure profiles on Figure 
12 indicate that the pressures increase when the gate seat is downstream 
from the design seat, and decrease when the gate seat is upstream. However. 
for gate openings of 8.8 and 23.5 percent of the design head, the pressures 
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When the gate seats upstream from the crest, there is a sharp reduction in ~ 
the pressure at the crest, particularly at the small gate openings. | 
The theory that a gate seat downstream from the crest will give a down- 
ward direction to the flow passing beneath the gate seems to be confirmed by 
these tests. Although the increase in pressure was small, it was noticeable 

for all crests. 

The lowest pressures occurred on the two crests that most nearly corres- 
ponded with the datum shapes, Figures 8 and 9, and generally occurred for a 
gate opening of about 33 percent of the design head. For the other three 
crests, Figures 10, 11, and 12, the lowest pressures occurred at the smallest 
gate openings. The lowest pressure recorded was about 20 percent of the de- 
sign head below atmospheric. 


Gate Seals 


Radial gates are usually provided with side and bottom seals. In some 
installations the side seal bears against a steel plate embedded in the face of 
the pier; this provides a smooth surface for the contact and prevents exces- 
sive wear of the sealing strip. The bottom seal contacts a metal plate set 
flush with the spillway surface; a shallow depression about 3/32 by 3 inches 
in the metal sill plate acts as a seating surface. 

In 1952, a limited number of tests were performed in the Bureau of 
Reclamation Hydraulic Laboratory to determine the effectiveness of the 
gates seals being used.(8) Because of limitations in the test rig, sliding of 
the side seal on the wall plate corresponding to the gate opening or closing 
was not possible. The seal could be moved in a direction normal to the wall 
plate, to determine the force required to seat the seal. When the seal seated, 
leakages were too small to be conveniently measured, so visual observations 
of the leakage sufficed for the tests. 

Two tests were made with each seal. The first test was to determine the 
lateral thrust per lineal foot of seal. If the lateral thrust was within reasona- 
ble limits, the seal was tested in the hydraulic rig to determine its action un- 
der pressure. 

The lateral thrust tests were made with a wooden mock-up of the clamp 
plates of a radial gate, Figure 13. The thrust was determined by lifting the 
seal up and off the laboratory floor with a spring scale and steel bar placed 
under and perpendicular to the seal axis. The seal was lifted at a point mid- 
way between two pieces of paper placed 12 inches apart on the floor under the 
seal. When the paper could be slipped free, the thrust load on the floor was 
considered relieved and the spring scale load recorded. The frictional drag 
was computed from this thrust value. 


Belt-type seal 

One proposed side seal was made from fabric-reinforced rubber belting. 
A seal 8 inches wide and 15 feet long was assembled in the wooden mock-up 
for lateral thrust tests. The lateral thrust was determined for various seal 
positions. Zero seal position occurs when the space between the end of the 
gate and the wall plate is three-quarters of an inch. The tests showed that at 
the zero seal position, the lateral thrust was 16 pounds per lineal foot of seal 
This thrust was high; and with the addition of a waterload at a 40-foot head, 
the frictional drag would be 75 pounds per inch of seal. In view of the high 
drag, no hydraulic tests were performed. 
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(A) General view of the 3/8 x 8 inch x 15 foot long 
belting seal bolted to a 40-foot radius, 


(B) Close-up view of the above assembly. 


FIGURE 13 
RADIAL GATE SEAL TESTS 
Wooden Mock-up of Rubber Belting Side Seal without Brass Shoes 
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Belt-type seal with brass shoes 

It was decided to investigate the same seal with brass shoes bolted to the 
belting to reduce the frictional drag, Figure 14A. The lateral thrust at zero 
seal position for this seal was 158 pounds per lineal foot; but because of the 
smaller coefficient of friction with the brass shoe, the frictional resistance 
at a 40-foot head was 34 pounds per inch of seal. 

Hydraulic tests showed that when the seal was moved away from the wall 
plate, the seal shoe slid downstream and remained in flat contact with the wa 
plate. However, when the gate was moved nearer to the wall plate, the brass 
shoe did not slide upstream on the wall plate, and the rubber belting buckled. 
Since this action was independent of the reservoir head, the seal was con- 
sidered unacceptable. 


Molded-type belt seal 

As a means of reducing the lateral thrust of the belt type of seal, it was 
proposed that the belt seal be molded to an angular shape, Figure 14B. The 
lateral thrust of this seal was negligible at zero seal gap. The frictional dra 
on the wall plate was computed at about 20 pounds per lineal inch. 

Hydraulic tests indicated that the brass shoe did not slip upstream freely 
as the seal gap was decreased. The seal buckled when the gap was decreases 
1/4-inch, and the sealing lip raised slightly from the seat. It was concluded 
that the seal operation was only fair. 


Angle seals 

Two types of angle seals were investigated. One was an all-rubber seal, 
and the second was similar, except that it was fabric-reinforced, Figure 15. 
The lateral thrust with these seals would depend on how they were mounted 
on the gate sides. Because of this, tests were performed to determine at 
what head these seals would seal and unseal for various seal-to-side wall 
spacing. In addition, identical tests were made for the all-rubber seal with 
a brass bar on the sealing edge. 

For the all-rubber angle seal, without the brass bar, the water head re- 
quired for sealing at a 1/2-inch gap was 22 feet, while unsealing occurred at 
an 8-foot head, Figure 16. With the brass bar, these values were 23 and 7 
feet. The fabric-reinforced seal sealed at a 19-foot head and unsealed at a 
6-foot head. The computed frictional drag for these seals was about 37 
pounds per lineal inch without a brass bar and 6.5 pounds per lineal inch witl 
a brass bar. 

Although the angle seal with the brass bar appeared to have better sealing 
and unsealing characteristics, as well as a lower frictional drag, it allowed 
slightly greater leakage than the all-rubber seal and probably would be less 
desirable for this reason. 

Two seals were made that included rubber-to-metal sealing surfaces and 
a brass member to reduce frictional drag, Figure 17. The seal with the rub 
ber flaps (Figure 17A) was generally unsatisfactory, since the head required 
to obtain a water seal at normal gaps was higher than normally found for 
spillway gates. 

The second seal, shown on Figure 17B, incorporated a more stabilized 
sealing lip. At a 1/2-inch gap, this seal sealed at a 20-foot head and un- 
sealed at a 4-foot head, Figure 18. The frictional drag with this type seal at 
a 40-foot head and at zero gap is approximately 11 pounds per lineal inch. 


RADIAL GATES 57 


(A) View of the rubber belting seal with brass shoes. 


(B) View of the belt type, molded angle seai, with brass shoe, 
cut from a fabric reinforced angle seal shown on Figure Sa. 


FIGURE 14 
RADIAL GATE SEAL TESTS 
‘ Belt Type Side Seals with Metal Shoes 
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(A) View of the angle seal with fabric reinforcement. 


“r= 


(B) View of the all-rubber angle seal with a brass bar added 
to reduce frictional drag on wall plate. 


FIGURE 15 
RADIAL GATE SEAL TESTS 


Views of the Molded Angle Side Seal with Fabric 
Reinforcement and a Brass Bar 
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(A) View of the angle seal with a brass bar plus a rubber sealing flap. 


(B) View of the angle seal with rubber removed from the lower leg 
and a brass plate added. 


FIGURE 17 


RADIAL GATE SEAL TESTS 


Views of the Angle Side Seal with a Brass Bar plus 
a Rubber Sealing Flap and with a Brass Plate. 
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Conclusions 


Based on the results of these tests and on the computations for frictional 
drag, a side seal design should include the features listed below. A suggest 
design incorporating these features is shown in Figure 19. 


1. The sealing force should not be dependent on the water head. 

2. A rubber-to-metal seal is preferred to obtain the least water leakage. 
The portion of the waterload carried by the rubber seal should be low so as 
to reduce the frictional drag. 

3. The heel end of the seal should be mounted so that it will not contact 
the wall plate and increase the frictional drag. 

4. The major part of the seal waterload should be carried through a metz 
member to the wall plates to minimize the frictional drag of the seal. 


Radial Gate Bottom Seal Tests 


The hydraulic test rig used in the side seal tests was modified to simulat 
a bottom seal arrangement, Figure 20. The seal seat was made with a curv 
depression 3/32-inch deep and 3 inches long measured parallel to the long 
dimension of the seal to represent a damaged seating surface which might b 
encountered in the field. The angle of the seal seat relative to the seal was 
similar to that of a field gate. Two all-rubber seal samples were available 
for test; each one was rectangular in section. One was made of stock with 
Type A Shore durometer hardness of 69; the other was made of stock with a 
hardness of 38. 

Both seal samples were tested in the hydraulic rig at heads up to 195 fee 
of water. This test consisted of compressing the seal against the simulated 
gate seat and increasing the hydraulic head to 195 feet of water. Both seals 
appeared to seal satisfactorily, since there was no leakage at the depressio. 
in the gate seat. However, at heads of 55 feet or more, the soft rubber seal 
had a tendency to creep under the downstream clamp plate. 

The principal test conducted on these two seal samples was to determine 
the force required to compress the seal. These data were desired for gate 
design purposes to assure sealing of depressions in the seal seat. This in- 
formation was obtained on a Universal-type testing machine by determining 
the vertical force required to compress a 1-foot wide section of the seal 
against wet and dry seats. The wet seal data were obtained by submerging 
the simulated seal seat in a pan of water. A plot of the compression force, 
pounds per lineal foot of seal for wet and dry seat, is shown for both seal 
samples in Figure 20. 


Conclusions 


The rectangular-type radial gate bottom seal tested was satisfactory. T 
softer rubber material would probably form a more effective seal against 
surface irregularities in the seal seat. The tests showed that for surface 
irregularities approximately 3/32 inch deep, the soft material would form z 
seal at an applied force of 50 pounds per lineal foot of seal; the harder ma- 
terial would require a force of 140 pounds per lineal foot. 
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NOTES 

1. Approximately one-third of seal 
water load across the 13° di- 
mension is carried by the gate, 
the remainder of water load is 
carried to the wall plate. 

2.Brass pipe may be slightly curved 
during assembly to fit gate 
curvature. 

3.At higher heads,seal leg will 
tend to curve,and thus releave 
the frictional drag of rubber 
sealing lip on wall plate. 

4.Drawing is to full scale, 


---% Holes, 13" Centers 
7g x-strong brass pipe 
4 123" wall thickness 


_--$" Initial deflection of 
seal on wall plate. 


RADIAL GATE SEAL TESTS 
PROPOSED SIDE SEAL WITH A BRASS PIPE 
AND A RUBBER SEALING LIP 


| FIGURE 19 
415 G 
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_ Test gate leaf 


Upstream clamp plate-------~~ 


_-- Depression, % x 3~inches wide, at angle 
iron & This cut represents an irregular 
in seal seat (sill plate) 


--Compression of seal 
~- 100% compression=0.575 inches. 


ni 


T RUBBER-38 DUROMETER 
TYPE A SHORE 


Force required to seal depression 
inseol seot.~. 
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1 
~HARD RUBBER—69 DUROMETER 
TYPE A SHORE 


COMPRESSION OF SEAL ~- INCHES 


700 
COMPRESSION FORCE-POUNDS PER FOOT OF SEAL 
FORCE REQUIRED TO COMPRESS RADIAL GATE BOTTOM SEAL 
RADIAL GATE SEAL TESTS 
ARRANGEMENT FOR HYDRAULIC AND COMPRESSION TESTS OF BOTTOM SEAL 
FIGURE 20 
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BOUNDARY-SHEAR STRESS IN UNSTEADY 
TURBULENT PIPE FLOW 


| M. R. Carstens,1 A.M. ASCE and John E. Roller,2 J.M. ASCE 


| INTRODUCTION 

There are two reasons for the preponderance of experimental and theoreti- 
jal knowledge concerning steady flow in circular pipes as compared with un- 
teady flow. First, problems involving steady or quasi-steady flow in pipes 
rise repeatedly in the majority of engineering works. Second, experimental 
esearch of unsteady flow has been delayed until the development of reliable 
ynamic-pressure-measuring instruments. 

In this paper, a theoretical analysis of boundary-shear stress for unsteady 
urbulent flow in a pipe is given, followed by a report of experiments with un- 
teady flow in a smooth pipe. 

There is good reason to restrict the study to turbulent flow. In laminar 
low, the instantaneous mean velocity and instantaneous velocity profile (hence, 
ne boundary-shear stress) are dependent upon the history of the motivating 
orce or piezometric-pressure gradient. As a result of this historical de- 
endence, each case of unsteady laminar flow is unique. Conversely, in tur- 
ulent flow, the strong lateral diffusion of the turbulent eddies tends to elimi- 
ate the dependence of the velocity distribution upon the history of the motion. 
‘this hypothesis is justified, then the dependent flow variable, boundary- 
hear stress, is a function of instantaneous values of the flow variables, 
ean velocity and acceleration (or piezometric-pressure gradient). Thus, a 
reater degree of generalization of results appears to be possible for un- 
teady turbulent flow in a pipe than for laminar flow. 


Ote: Discussion open until July 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1945 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 2, February, 1959. 
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NOMENCLATURE 


cross-sectional area of pipe; 


2 
1 f (>) dA, linear-momentum-correction coefficient; 
oa 


inside diameter of pipe; 

Darcy-Weisbach coefficient or boundary-drag coefficient; 
acceleration of gravity; 

piezometric head; 

piezometric head in the reservoir; 

length of pipe 

reciprocal of the exponent in the power law of velocity distribution; 
piezometric pressure; 

discharge or volume rate of flow; 

radial coordinate measured from the centerline of the pipe; 
pipe radius; ) 

fluid density; 

shear stress; 

boundary-shear stress; 

time, measured from the instant of initiation of the flow; 
velocity in the pipe; 

mean velocity in the pipe; 

“observed” velocity determined from the jet coordinates; 
volume of a fluid element within the pipe; 

axial coordinate measured from the upstream end of the pipe; 
subscript designating unsteady flow; and 


subscript designating steady flow. 


Theoretical Analysis 


The linear-momentum equations are the beginning point of a rational ana 


sis of boundary-shear stress. In uniform unsteady flow in a pipe, the veloci 
v, is everywhere axial and is a function of r and t only. Similarly, the she 
stress, T , is a function only of r and t. Since there is no mean transverse 
acceleration, the piezometric pressure, PS , is a function only of x and t. 


The fluid element, of length dx, extends to the pipe wall, as shown in 


Figure 1. The linear-momentum equation for this element is 


— opr. dx Bear aon ot ee gee joes + 2 (evds 
ox Ox iN ot aa 
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T 27 bs dx 
fo AY et 2 ES 1 


JK» + 


Figure 1. Fluid Element. 


The terms on the left side of the equation are the external forces on the 
fluid element. The first term is the combined pressure and weight force in 
the x direction. The second term is the boundary-shear force. 

The first term on the right side of the equation is the difference of momen- 
tum flux through the end areas of the fluid element at a given time. Since the 
equation is being applied to uniform flow, this term is zero. The existence 
of turbulent velocity fluctuations has no effect on the magnitude of this term. 

The second term on the right side of the equation is the time rate of change 
of linear momentum of the particles within the interior of the element. Since 
dx is independent of r and since v is not a function of x, the following sim- 
plifications are possible, 


= = erent dv dx 
) d¥ Sach ie vant dx = edQ dx = erg tT. ah 
2 fev ot - dt dt 


Avain the existence of turbulent velocity fluctuations has no effect on the 
nagnitude of the term. 

_ The simplified ites Mnomentin equation for the fluid element of Figure 1 

S 


* 
: eee aEey 
om a dt (1) 


_ Another linear-momentum equation can be derived in a similar manner 
vith an annular cylindrical element of radius r and thickness dr. The result- 
ng linear-momentum equation is 


1 
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=O) nes = fo) 
Opy es otra) cae (2) 


Ox Ve 


In Equation (2) Op*/Ox , is not a function of r; but, since v is a function of 


r, Ov/Ot isa function of r. Therefore the second term on the left side of 
Equation (2) is a function of the velocity distribution. In contrast, the cor- 
responding term for steady flow is a constant regardless of the velocity 
distribution. 

If the magnitude of the acceleration is small, a reasonable assumption is 
that the velocity distribution is the same for the unsteady flow as for the cor- 
responding steady flow. The basis for the assumption is the strong lateral 
diffusive effect of the turbulent eddies. The process can be visualized by 
considering the chronological sequence of events as an initially steady flow is 
accelerated to another steady-flow state. Any deviations from the steady- 
state velocity distribution which occur during acceleration will be rapidly 
removed by the turbulent-eddy diffusion in the final steady state. In addition, 
the turbulent-eddy diffusion is also present during the acceleration. Thus, 
for turbulent flows undergoing moderate acceleration, equality of the unsteady 
state and steady-state velocity distributions is a rational assumption. 

The power law of velocity distribution for turbulent flow in a smooth pipe 
is used in the subsequent analysis. The power law of velocity distribution is 


(1). 


vy = (2n¢)(nel) (1 = 1 WO 
V 2 ine ) 


0 


The value of n is seven for pipe Reynolds numbers less than 1 (105). The 
value of n increases with increasing Reynolds number, becoming 10 ata 
Reynolds number of 3.2(106). 

Equation (2) is integrable with respect to r upon substitution of the power 
law of velocity distribution. The piezometric-pressure gradient is eliminates 
by means of Equation (1). The resulting expression for shear-stress is 


hte a Tou(=) t [F(t] e dv (3) 
and for shear-stress gradient is 
or = emt Fe(=) egy (4) 
r Fo % t 


Values of Fy (r/rp) and Fo9(r/ro) are given in the Table for n= 7. The 
shear-stress gradient is a function of r in contrast to steady flow in a pipe. 
The unsteady-flow shear-stress in Equations (3) and (4) would be more 
meaningful if related to the well-established expressions for shear stress in 

steady flow. This relation is possible by virtue of the equality of the shear 

stress and the shear-stress gradient at the centerline for both steady and un. 
steady flow. Since similar velocity distributions were assumed in steady anc 
unsteady flow, the deviation of unsteady-flow shear-stress distribution from 
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inear must be solely the result of acceleration. If all fluid particles were 
jaccelerated equally, as in rigid-body motion, the internal stresses (shear 
stresses) would not be influenced by the acceleration. Therefore, the change 
in the shear-stress distribution must be the result of relative acceleration 
between adjacent particles in the cross section. Inasmuch as the velocity 
istribution is symmetrical about the pipe axis and is continuous through the 
axis, the relative acceleration between a particle on the axis and a particle 
an infinitesimal distance from the axis will also be infinitesimally small. 
Thus, as a consequence of the assumption of identical velocity distributions 


in unsteady and steady flow, the shear-stress gradients at the axis are equal, 


r=0 Fr Ire0 


In addition, the shear stress is zero at the pipe axis for either steady or un- 
steady flow, or 


T. si dF = © 
Ss 0 ee 0 
But in steady flow 
aia Sa elos 
dr ae 
Introducing these concepts into Equation (4), 
z v3 dV 
toads Fair \ . ze (5) 
9 ro Pol ir =O 
Similarly, for Equation (3), 
Ty = Tos(Z | + [Fa(e\|ro e dV (6) 
es rs dt 


Values of F3/r/r) are given in the Table for n= 7. The physical concepts 
embodied in Equations (4), (5), and (6) are more easily visualized by referring 
to Figure 2. 

Equation (5) is the desired relationship between the boundary-shear stress 
in unsteady flow and that in steady flow. Rearranging Equation (5) and intro- 
ducing the boundary-drag coefficient, f = 8T)/P v2, results in 


= 1+ 0449 _D,- dV 7 
fu aS Viewed t ey 


The boundary-drag coefficient, fy, is influenced in two ways by the magni- 
ude of the Reynolds number. First, fg is a function only of the Reynolds 
4umber for flow in a smooth pipe. Second, the value of n, which characterizes 


4 
; 


. 
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ou (accelerated) 
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y (decelerated) 


oe (accelerated) 


an (decelerated) 
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Figure 2. Shear Relationships in Unsteady Turbulent Flow. 


the shape of the velocity profile, is also a function of the Reynolds number. 
The numerical value, 0.449, in Equation (7) is for n = 7. This numerical vah 
is 0.391, 0.346, and 0.310 for n of eight, nine, and ten, respectively. 

Thus, the theoretical analysis is indicative that the boundary-shear para- 
meter, f,, of unsteady and turbulent flow in a smooth pipe is independent of 
the history of the motion and is dependent only upon the corresponding stead: 
state value of fg and the instantaneous value of the acceleration parameter, 
(D) (dv/dt) / V2. Equation (7), which is a mathematical expression for this_ 
dependence, was formulated using the hypothesis that the velocity distributio 
is similar for both steady and unsteady flow. For moderate acceleration, th 
hypothesis appears to be valid, but experimental results are desirable in or« 
to verify and to establish the limits of the hypothesis. ; 
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TABLE: Values cf F F,, and F LOweN Tae) | 


eels a 3 
r/r F, (r/r,) F,(r/r,) F(r/r,) 
| Equation {3) Equation (4) Equation (6) 
| 0) 0) -0.1123 ) 
Oe =~ 01: -0.0962 0.0002 
0.2 -0.0200 -0.0861 0.0025 
0.3 -0.0279 -0.0706 0.0058 
0.4 -0.0340 -0.0534 0.0109 
O25 -0.0385 -0.0320 0.0176 
0.6 -0.0403 -0.0080 0.0261 
Ont -0.0396 0.0256 0.0378 
0.8 -0.0349 0.0706 0.0549 
0.9 -0.022} 0.1414 0.0786 
LO ) 1.0000 0.1123 


Experimental Study 


Whether steady or unsteady flow, experimental determination of boundary- 
shear stress is indirect. In unsteady flow, this determination, by means of 
Equation (1), is complicated by the necessity of measuring p*(x,t) and V(t). 
in these experiments, the function V(t) was determined by analysis of a 
motion-picture record of the free jet of the pipe outlet. Commercially avail- 
able pressure transducers were used to determine p*(x,t). 

The boundary and initial conditions of the experiment were as follows: 

a) the smooth, straight, circular pipe was attached to a large reservoir of 
water in which the piezometric head was maintained constant; and (b) the 
low was initiated by the instantaneous opening of a downstream valve. Since 
he fluid was at rest prior to the opening of the valve, the initial flow was 
aminar throughout the pipe. With the passage of time, spots of turbulence 
vere generated. These turbulent spots increased in size with increasing 
ime until the flow was turbulent throughout the pipe. 


ixperimental Equipment and Technique 


The overall arrangement of the equipment is shown in Figure 3. The test 
ection was extruded 1/2-inch brass pipe. The experimentally determined 
iezometric-head gradient of six steady-flow runs was within the range, +0.1 
er cent to -1.08 per cent, of the piezometric-head gradient computed from 
he Blasius equation. This close agreement means that the pipe was smooth 
nd that the diameter, measured by means of a micrometer, was correct. 
iezometers were spaced at intervals of 95 diameters along the pipe. The 
ressure-transducer signals were amplified electronically and were recorded 
n a Sanborn two-channel oscillograph. The configuration of the downstream 
st from the horizontal pipe was recorded by means of a 35 mm motion- 
icture camera. i 
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erimental Data 


The oscillographic records of pressure were replotted as instantaneous 
siezometric-head lines in order to delineate the laminar - and turbulent-flow 
egions. The succession of piezometric-head lines of a typical run are shown 
in Figure 4. The abscissa of Figure 4 is the distance along the pipe axis 
rom the reservoir. The ordinate is h/hg. The position of the zero value of 
€ ordinate is lowered with increasing values of tV gh /\ in order to 


void the superposition of many piezometric-head lines. The dashed line in 

igure 4 delineates the regions of laminar and turbulent flow. Laminar-flow 

egions are above the dashed line and turbulent-flow regions are below the 
hed line. 

With the free-fall equations, the “observed” velocity-time relationship of 
Ka downstream jet was determined from the coordinates and from the time 
recorded on the 35 mm film strip. Two assumptions are implicit in the use 
of the free-fall equations. First, air resistance is negligible on the jet. 
Second, the jet is composed of individual masses in which each fluid particle 
has the same linear momentum at the point of measurement. Two experimen- 
tal observations substantiated the validity of the assumptions, as follows: 

(a) the “observed” velocity-time relationships were identical whether the 
point of measurement was 0.5 ft below the pipe outlet or whether the point of 
measurement was 2.50 ft below the pipe outlet; and (b) in all the film strips, 
discrete masses of fluid appeared to be forming in the jet at the point of 
measurement. However, the “observed” velocity, V', is not the mean veloci- 
ty in the pipe. 

The mean velocity, V, in the pipe is less than the “observed” velocity, V', 
determined from the jet coordinates. The difference in the two velocities is 
due to the radial variation of velocity within the pipe. The radial variation of 
velocity at the pipe outlet is the result of the boundary-shear force as the 
fluid passes through the pipe. As the fluid passes the pipe outlet, the boundary- 
shear is decreased suddenly to a negligible value and the internal shear 
forces within the jet tend to eliminate the non-uniform velocity distribution. 
The linear momentum of the jet is greater at the pipe outlet than if the veloci- 
ty distribution were uniform. The mass of fluid passing the pipe outlet ina 
time at is PQat. The linear momentum of this fluid mass is Cm PQVat. 
The linear-momentum-correction coefficient, Cy,, is a measure of the non- 
iniformity of the velocity distribution as the fluid passes the pipe outlet. In 
he absence of any horizontal external forces the linear momentum of this 
mass is constant in the jet. At the point of measurement of jet coordinates, 
he linear momentum is PQV' at. Therefore, the desired quantity V is 
squal to V'/ Cy: Only the value of the “observed” velocity, V', can be de- 
ermined from the film strip of the jet. Other methods had to be employed to 
letermine Cy. 

The experimentally determined values of Oh/Ox of laminar flow were in- 
-orporated into the solution of the Navier-Stokes equations in order to deter- 
nine the velocity distribution at the pipe outlet as a function of time. It is to 
e noted in Figure 4 that the flow was laminar at the pipe outlet for 
Vgho/L <1.38 even though turbulent flow existed in portions of the pipe 
ta much earlier time. The heat diffusion equations analogous to the Navier - 
ttokes equations have been solved with analogous boundary and initial condi- 
ions and with the analogous piezometric-head gradient as an arbitrary func- 
ion of time.(2) The piezometric-head gradients were determined from 
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Figure 4. Piezometric-Head Lines of a Typical Run 
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igure 4 and were replotted in Figure 5. A simple analytical function, 
h/Ox ~ cos Kt was selected as being a good approximation for the experi- 
entally determined piezometric-head gradient of laminar flow. Of course, 
- had to be determined empirically for each of the eight runs. This function 
as then incorporated into the solution of the Navier-Stokes equations for the 
ermination of C,, as a function of time. Because the length of the 
inar-flow region was less than the distance between pressure measure- 
ent stations when tVgho/L >1.00, the piezometric-head gradient of - 
inar flow could not be experimentally determined from Figure 4 when 
Vgho/L >1.00. Thus during the interval 1.00 < t-/gho/L < 1.38, Cm had 
be determined in some other manner. 
The magnitude of the discontinuity of V' at tVgho/L = 1.38 shown in 
igure 6 was utilized to obtain the value of Cy, at this time. The discon- 
inuity occurs at the time the turbulent flow first passes the outlet. The dif- 
erence in velocity distribution between the turbulent flow and preceding 
aminar flow results in a discontinuity of linear momentum in the laminar 
und turbulent jets at this time even though V is continuous. The value of 
fh of the turbulent flow was assumed to be 1.02 based upon the seventh-root 
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Figure 5. Piezometric-Head Gradients of a Typical Run. 
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law of velocity distribution. Since V' must equal C,,V and since V is 
equal for both laminar and turbulent flow at this time, the experimental 
point, shown as a circle in Figure 6 was readily determined. The dashed 
line in Figure 6 is an interpolation between the two methods. 

The resulting V-t function is shown in Figure 6.* 


Experimental Results 


The boundary-shear stress was determined indirectly by means of 
Equation (1) utilizing the data shown in Figures 5 and 6. The results are 
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Figure 6. "Observed" Velocity, Mean Velocity, and Linear -Momentum 
Coefficient of a Typical Run. : 


*A more complete discussion of the method of analyzing the data is containe 
_in the paper “Transition from Laminar to Turbulent Flow in a Pipe,” by 


M.R. Carstens, Journal of the Hydraulics Division, Proc. A.S.C.E., Paper 
1450, December, 1957. 
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fesented in Figure 7 in the form suggested by Equation.(7). Results of ex- 
iments conducted at Massachusetts Institute of Technology(3) in an 
—o water tunnel are also shown in Figure 7. 


salysis of Results 


| There are two obvious explanations for the scatter of the experimental 
sults shown in Figure 7. 

| First, the scatter of the experimentally determined values is an indication 
the precision of determining the drag coefficient, f,. The p-t function 

uld be measured within a probable error of + 2 per cent in regard to pres- 
re. In view of the complexities in determining the V-t function from the 
tion-picture data, a probable error of + 3 per cent would seem to be 


aia 


a 


dary-Shear 
Fi e 7. Experimentally Determined Values of Boun 
Ps Stress of Unsteady Turbulent Flow in a Smooth Pipe. 
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optimistic. A further source of possible error is in the graphical determin 
tion of dV/dt from graphs similar to Figure 6. Since none of these pessib 
errors is compensating, the entire scatter of the experimental results migh 
be attributed logically to experimental error. 

Second, the unsteady-flow drag coefficient, f,, may be a function of the 
history of the motion. The presentation of f,/fg as a function of 
D(dV/dt)/f,V2 is deduced from Equation (7). However, the basic assumptio 
in the derivation of Equation (7) is that the history of the turbulent flow does 
not influence the instantaneous flow variables. Thus, if history of the unstea 
motion is important, the scatter of the experimental results in Figure 7 mig 
be the result of omitting the parameters necessary to define this history. Ir 
this regard, Daily, et al(3) noted a historical effect during deceleration but 
not during acceleration. 


CONCLUSIONS 


The experimental results are not sufficiently precise to either prove or 
disprove the validity of the assumption that the history of turbulent flow doe: 
not influence the instantaneous flow variables. The derivation, leading to 
Equation (7) for the drag coefficient of unsteady flow is considered to be ra- 
tional. 

Both the results of this study and those of Daily, et al,(3) indicate that the 
current engineering practice of assuming that f, =f, is acceptable for turk 
lent flow in pipes. In the first place, the relative importance of the boundar 
drag force decreases in unsteady flow as compared to steady flow while the 
disturbing force and the inertial reaction increase in relative importance. ] 
the second place, the unknown resistance and virtual-mass effects of variou 
fittings in the pipe line are likely to be of greater importance than the small 
changes in f, indicated by Equation (7). 
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| CHARACTERISTICS OF FLOW OVER TERMINAL WEIRS AND SILLS? 


Corrections to Discussion by Corrado Ruggiero! 


CORRECTIONS~—In printing F. Paderi’s discussion of this paper in 
iroceedings Paper 1832 (October, 1958) the fifteen references used in that 
iscussion were inadvertently omitted. We offer them at this time for the 
enefit of interested readers: 


All these branches of research offer several developing possibilities, as 
- is attested by the authors surveys, which are very valuable on either the 
xperimental side or because of reflections on the theoretical field. 
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THE HYDRAULIC JUMP AT AN ABRUPT DROP@ 


Closure by Walter L. Moore and Carl W. Morgan 


| 
| 
| 
| 


WALTER L. MOORE} A.M. ASCE and CARL W. MORGAN,2 A.M. ASCE.— 
The authors wish to thank Mr. McPherson and Mr. Forster for their valuable 
discussions. 
_ Mr. McPherson has compared and integrated the findings with those of 
others. He asks about the influence of the direction of movement of the tail- 
water level on the limiting depths for the various types of jumps. A brief 
investigation of the effect of rising and falling tailwater, made at the time of 
the first experiments, indicated a negligible effect of this variable provided 
sufficient time was allowed for equilibrium to be established. As a result of 
the comments by Mr. McPherson and the data he presented additional ex- 
periments on this effect were made in the Hydraulic Laboratory at the Uni- 
versity of Texas. These experiments indicated that if sufficient time is al- 
lowed the effect of the direction of tailwater movement was of the same 
magnitude as the experimental error in determining the tailwater depth for 
Co) 
Y4 
= 4 and Fy =6 it was possible to set the tailwater so that the 


Froude numbers less than 6. In fact for = 3 and Fy =4 and also 


for 


fe) 
¥y 
jump would slowly alternate between the B type and the wave type. In order 
to change from the B type to the wave type or vice-versa in a moderate time 
it was necessary to set a higher value with rising tailwater than with falling 
tailwater, a result similar to those observed by others. Those results, how- 
ever, were for a different range of Fj and relative drop height AZ,/Y¥4 than 
used in the authors’ experiments. The equilibrium condition which the 
authors obtained for F4<6 could not be obtained at Fy = 8. For AZ,/Y1 =4.0 
za a rising tailwater the change from the B type to the wave type required a 
a = 12.3 while for a falling tailwater the change occurred at = = 11.4. 
it is felt that this may have been due to the fact that at this Froude number the 
op of the wave rose higher than the walls of the flume. On the other hand it 
may be that at this high Froude number the flow pattern is such that the 
‘hysteresis effect” mentioned by McPherson and others is present. 

In the additional experiments mentioned above, tailwater depths were as- 
sertained by means of piezometers rather than with a point gage as in the 


: , (0) 
riginal experiments. The experiments wi Bee = 3 were made at a scale 


1. Proc. Paper 1449, December, 1957, by Walter L. Moore and Carl W. 
Morgan. 

.. Chmn., Dept. of Civ. Eng., Univ. of Texas, Austin, Tex. 

}. Asst. Prof., Dept. of Civ. Eng., Univ of Texas, Austin, Tex. 
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twice that used in the original experiments. These new tests indicated that 
the values shown in Fig. 7 for the maximum jump B are slightly high. Fig. Ta 
is here presented to better define the several values required to insure a 
jump B or a wave to form. 

Mr. McPherson calls attention to the disagreement with data presented by 
Hsu(4,5) in which a systematic transition between the types of jump (A,B, 
Wave) was found as a function of the Froude No. Fj. Actually the investiga- 
tion reported by the authors was started to explore some longitudinal charac- 
teristics of the jump at a drop assuming the existence of the systematic 
transition reported by Hsu. It was soon found that type B jumps could be ob- 
tained at low values of Fj where Hsu’s curves indicated only jump A would 
form. Further investigation showed that the full range of jump types was 
possible at any value of Fy (at least above 2.0) controlled only by the tail- 
water elevation. The systematic transition as a function of Froude No. was 
not observed by the authors. 

From considerations similar to those leading to figure 13 it is apparent 
that a tailwater control which produced a relatively flat curve of tailwater 
elevation vs discharge as opposed to the steep one shown in Fig. 13 would 
result in a systematic transition from Jump B through the wave to Jump A 
as the discharge was decreased. With the depth Yj held constant a reduction 
in discharge would correspond to a reduction in Fy. This could explain a 
systematic change of jump type with decreasing Fj as reported by Het, 
However, it is not clear how the slope of the curve of tailwater depth a vs 
Froude No. Fy could be positive for the type A & B jumps and negative for 
the wave type jump as indicated in his Fig. 10 of Reference 5. 

Mc. McPherson’s Fig. D-2 calls attention again to the fact that the band 
shown for Jump A in Figs. 6, 7, and 8 does not represent the maximum 
values of the tailwater. In the Results section of the paper it is stated that 
“The conditions under which Jump A will form are shown as a very narrow 
band, but actually Jump A will form at values of Yo/Y¥4 greater than the 
upper bound------ ce 

In response to Mr. McPherson’s suggestion values of Yo/Y4 are listed 
below for the various forms of the jump for which relative bottom velocities 
are presented in Figs. 9, 10, 11. 


Mr. Forster effectively emphasizes the importance of the type of jump 
(type A, wave, or type B) in determining the severity of attack on the channel 
bottom just downstream from the drop. His examples well illustrate the ef- 
fectiveness of Jump A and the wave in protecting the bottom from high veloci- 
ty currents. However, the wave type can result in some fairly high bottom 
velocities for large drops ee = 3; ae = 4 Figs. 10 and11). In 

: 1 
Figures 9, 10, and 11 the various regions and curves are identified by a letter 
and a number as listed in the table. It can be seen that with the wave type and 


a ales number of 4 a downstream velocity of 3V59 may be attained for 
te) AZ . 
Lee = 3.0 and 4V> for = = 4.0. It may seem odd that a Froude number 


of 4 resulted in higher relative bottom velocities than Froude numbers of 2, 
6, and 8. This does not seem unreasonable, however, when the complex flow 


pattern below the drop is considered. With a large 


and a small Froude 
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Z 
— -2 Fig.9 
ay 
YQ %2 Y2 
sabe — eF, — 
Type F; Y Type Fy % Type Fy V1 
Bi2meoak B4 5.4 605-4 
W2 4.3 Ww4 6.8 Wé6 9.4 
A2 4.6 A4 7.1 A699 
Z 
) : 
en ee : 
Yi Fig. 10 
B2 4.1 B4 6.3 B6 8.4 B8 11.4 
W2 5.1 W4 7.2 W 6 10.4 W 8 12.4 
A2 5.4 A4 17.8 A6 10.6 A 8 13.5 
oO : 
— =4 Fig. ll 
Yi 6 
B2 4.2 B4 5.9 B6 8.80 B 8 10.5 
W2 6.3 W4 8.5 W 6 10.7 W 8 11.8 
A2 6.5 A4 8.9 A6 11.4 


number the momentum of the jet simply is not sufficient to penetrate down to 
the bottom. At an intermediate Froude number the jet falling down from the 
wave still has sufficient momentum to penetrate down to the bottom. At 
higher Froude numbers the rise in tailwater depth is sufficient to lift the jet 
off the bottom, thus preventing high velocities from developing on the bottom. 
Mr. Forster points out that an overdimensioned drop may be unconserva- 
tive in that Jump B may be produced with associated high bottom velocities. 
This is certainly true for the conditions he assumed; that is, that Y9/Y, is 
constant. This is equivalent to assuming that the size of the drop is con- 
trolled by the elevation of the upstream channel while the downstream channel 
elevation and tailwater depth remain fixed. If, on the other hand, the size of 
the drop were controlled by lowering the bottom elevation of the downstream 
channel for a suitable distance while the tailwater elevation remains constant 
as controlled by conditions farther downstream, an oversize drop will then be 
conservative. Using the same initial conditions as in Mr. Forster’s example 
(Fy = 2.8 and Yo/Y, = 5) and a drop of AZ, = 2Y 1, Fig. 6 will give a minimur 
wave type jump (bordering Jump B). If the oe is increased to AZ, = 3Y}j, by 
lowering a suitable length of the downstream channel bottom, the values be- 
come Fy, = 2.8, Y2/Y1 = 6, yielding from Fig. 7a a maximum wave type jump 


(bordering on Jump A). In this manner the low bottom velocities would pre-— 
served even with an oversize drop. 
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Although Figs. 9, 10, and 11 give some very valuable data indicating the 
range of velocities to be expected near the channel bottom, they are not in- 
tended to be used for precise estimates of velocities in particular instances. 
Due to the complex nature of the flow pattern there is necessarily considera- 

| ble scatter in the data. In addition the scour potential of the flow is certainly 
dependent upon the turbulence characteristics of the flow as well as the mean 
velocity. Thus the velocities should be thought of merely as an index of the 
| Beour potential for comparisons between the various types of jumps. 
. Mr. Forster points out in his example that Jump B produces a bottom 
Felocity of the same order of magnitude as a drowned jump on a flat apron. 
|The peak values of the bottom velocities in Figs. 9, 10, and 11 are of the 
order of 6V> which is the same as that given by Mr. Henry (Trans. ASCE, 
1950, p. 690, Fig. 33) for a moderately drowned jump. These relative high 
velocities are not necessarily harmful if a paved apron is provided down- 
stream from the drop for a sufficient distance to protect the channel bottom. 
For a situation in which a rising discharge and consequent rise in tailwater 
would cause increasing submergence of the jump as its size increases, the 
introduction of an abrupt drop would be beneficial by preventing severe drown - 
ing. Instead of drowning the jump would shift toward the wave and Jump A 
with consequent reduction of attack on the channel bottom. 
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FLOOD FREQUENCIES DERIVED FROM RAINFALL DATA2 


Closure by J. L. H. Paulhus and J. F. Miller 


(1958) Journal of the Hydraulics Division, in Proceedings Paper 1856, there 
appears a discussion by C. O. Clark of Proceedings Paper 1451 that the 
authors were unaware of when they prepared their closing remarks to the 
discussion by Messrs. Riggs and Benson (Proceedings Paper 1690, June, 1958). 
The following remarks, therefore, are the closing remarks to Mr. Clark’s 
discussion and form an addendum to the closure published in November: 

Generally speaking, the differences between Clark’s opinions and those of 
the authors are of a philosophical as well as statistical nature. The writer 
criticizes the authors’ use of the term “unusual” in describing the 1933 
Codorus Creek flood. Clark claims that the most intense 6-hour rainfall 
over 10 square miles in the storm producing this flood, i.e., 6.8 inches, is 
not a sufficient approximation to the estimated probable maximum value of 
24 inches to be called unusual. Is that a measure of unusualness? 

Methuselah is reputed to have lived 900 years. Hence, following Clark’s 
reasoning, it would not be unusual for someone to live 200 years! 

The results of a recent extensive rainfall-frequency study(1) involving all 
available precipitation records in the Middle Atlantic region indicate that the 
shortest average recurrence interval that could reasonably be assigned the 
6.0 inch, 10-square mile, 6-hour rainfall of the October 1942 storm for either 
Codorus or Passage Creek basins would be about 160 years. Now, there is 
no reason why a rainfall of that intensity should always occur when the sur- 
face layers of the soil are, practically speaking, saturated, as they were in 
the case of the October 1942 Passage Creek flood. The same rainfall inten- 
sity could occur under much less favorable flood-producing soil conditions. 
Consequently, the shortest average recurrence interval of the October 1942 
flood can reasonably be expected to be longer than 160 years, indicating that 
the writer’s flood-frequency curve for Passage Creek probably yields ex- 
cessive peak discharges for the longer recurrence intervals. Examination of 
the streamflow records reveals that the curve very likely yields excessive 
discharges for the shorter recurrence intervals also. For example, the curve 
shows a two-year peak discharge of about 6000 cfs although that value has 
been equaled or exceeded only three times in the 26-year period, 1932-57 (see 


Basin D in Table 1 below). 


| 
ig ls; H. PAULHUS,! M. ASCE, and J. F. MILLER2.—In the November 


a. Proc. Paper 1451, December, 1957, by J. L. H. Paulhus and J. F. Miller. 
1. Staff Hydrologist, U. S. Weather Bureau, Washington, D. C. 
2. Meteorologist, U. S. Weather Bureau, Washington, D. C. 
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The shortest average recurrence interval that might reasonably be as- 
signed the 7.4 inch, 10-square mile, 6-hour rainfall of the 1899 storm men- 
tioned by the writer is about 800 years. The peak discharge of 36,000 cfs 
that Clark estimates could happen on Codorus Creek from such a rainfall 
under optimum flood-producing conditions would very likely have a longer 
average recurrence interval than 800 years. It should be realized, however, 
that even if the average return period of such a flood were 1000 years, 1 per- 
cent of the intervals between occurrences might be less than 10 years, and 5 
percent might be less than 51 years.(2) Obviously, the occurrence on 
Codorus Creek of a peak discharge of this magnitude within the next few 
years should not be astonishing. Nevertheless, the event may still be con- 
sidered unusual. 

Now, how about the unusualness of the August 1933 maximum peak dis- 
charge of 11,200 cfs observed on Codorus Creek? This was still the maxi- 
mum observed for the 25-year period ending with water year 1957. Accordin 
to the laws of probability there is a 1-percent chance that the average return 
period of this flood peak could be over 2000 years, and a 5-percent chance 
that it could be about 500 years.(2) Consequently, the 11,200 cfs peak dis- 
charge could be an unusual event. 

The chief trouble with Clark’s line of reasoning is that there is no objec- 
tive way for determining when or where to stop. For example, the Smethport 
Penn., storm of July 17-18, 1942 (OR 9-23) produced the maximum 10-square 
mile, 6-hour rainfall ever observed in the United States, namely, 24.7 inches. 
For some reason, the writer does not mention this storm. Perhaps it was 
believed that rainfall of that intensity could not occur over either Codorus 
Creek or Passage Creek basins. Why not? The probable maximum rainfall 
for both basins is estimated to be somewhat higher than for the storm site. 
Presumably then, rainfall of equal or greater intensity than for Smethport 
storm could occur or could have occurred over these two basins. The fact 
remains, however, that there is no evidence of the occurrence of such a 
rainfall on either of these basins within the period of our recorded precipita- 
tion or streamflow history, and the authors know of no way to base frequency 
analysis on what “could happen.” 

Referring to Passage Creek, Clark states, “This stream had its floods; 
Codorus Creek did not. At least not in the period 1933-54. In the 26-year 
period 1932-57, Passage Creek experienced three peak discharges exceeding 
the maximum for Codorus Creek. Is it not possible that Codorus Creek is 
less prone to flooding than Passage Creek? On the other hand, is it not pos- 
sible that Passage Creek, through chance distribution, has received more 
than its share of floods for this 26-year period? Table 1 suggests that the 
latter is definitely possible. Passage Creek is the only basin for which the 
three highest annual flood peaks are 3 to 5 times the fourth highest. A fre- 
quency analysis (Gumbel method adjusted for partial duration series) of its 
annual floods for the 26-year period, 1932-57 yields the following: 


Return period (years) 2 5 10 25 50 100 
Peak discharge (1000 cfs) 4.5 8.5 11.5 15.5 18.6 21.5 


These data differ from Clark’s frequency curve, presumably based on re- 
gional considerations, by less than 10 percent except for return periods 
shorter than about 5 years. Comparison of the statistical arrays of annual | 
floods for the seven basins of Table 1 indicates that the Passage Creek recor 


2 DISCUSSION 


Table 1. STATISTICAL ARRAY OF ANNUAL PEAK DISCHARGES (100 CFS) 


FOR SEVEN BASINS OF 61 to 101 SQUARE MILES 


Rank Basin 
A B & D E F G 
a eee ee 
1 97 112 41 210 150 160 97 
2 9 61 41 147 3 137 92 
5) 50 re) 36 as 65 115 19 
4 43 he 34 41 36 112 TT 
5 he 32 34 39 29 110 73 
6 ko 32 30 38 28 110 69 
7 4a 1. 30 37 27 95 66 
8 38 27 30 32 26 70 59 
9 35 27 29 31 26 70 oT 
10 33 25 28 30 2k 70 47 
11 29 2h 2T 26 23 68 39 
12 27 22 25 26 23 64 35 
13 26 22 a4 23 23 64 32 
14 2k 19 22 21 22 62 32 
15 2k 19 21 20 22 62 26 
16 24 18 20 18 21 60 18 
aly ¢ 22 18 19 18 20 kg 17 
18 21 17 18 17 20 39 17 
19 21 16 17 16 20 38 
20 19 16 16 14 ay 36 
21 17 14 16 14 17 32 
22 16 14 16 13 15 27 
23 15 11 11 12 14 
24 14 11 12 12 
25 14 9 6 We 
26 8 ) de 
27 10 


Chester Creek near Chester, Penn., 61 sq. mi., 1932-57 
Codorus Creek at Spring Grove, Penn., 74 sq. mi., 1933-57 
Linganore Creek near Frederick, Md., 82 sq. mi., 1935-57 
Passage Creek at Buckton, Va., 87 sq. mi., 1932-57 
Seneca Creek at Dawsonville, Md., 101 sq. mi., 1931-57 
Tohicken Creek near Pipersville, Penn., 97 sq. mi., 1936-57 
Tye River near Lovingston, Va., 92 sq. mi., 1939-57 


° 
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provides about the most unrepresentative sampling for the group. Can : 
Clark’s curve then be accepted as representative for this group of basins? 

The authors were unable to determine the plotting formula used or the 
criteria for construction of the curve in the writer’s frequency plot. Neither 
could it be determined why the frequency curve ascribed to the authors was 
plotted from their Table 4. For the Codorus Creek flood peaks plotted, the 
curve should have been drawn from data in the authors’ Table 3. Neither 
set of data is recommended as representative for Codorus Creek. For the 
situation under discussion, the authors recommend the use of all available 
streamflow records for the basin supplemented as much as possible by syn- 
thesized flood peaks. However, the synthesized flood peaks for the full 
period of rainfall record should yield good results. Table 2 compares these 
results with those from the observed record and shows that the difference is 
appreciable. The data were taken from Tables 3 and 6 of the authors’ 
original paper. 


Table 2. FLOOD FREQUENCIES (100 CFS) — CODORUS CREEK 
(By Gumbel Analysis — Not Adjusted for Partial Duration Series) 


Return Period (years) 


2 5 10 25 50 100 
Observed record, 1933-54 25 48 64 83 97 112 
Synthetic record, 1900-54 22 38 49 62 72 82 


In view of Clark’s questioning the reliability of flood frequencies derived 
from something other than floods and from data outside the watershed, the 
authors were astonished by his reference to rainfall data outside the subject 
basins as an indication of the possible frequency of extreme values of ob- 
served peak discharges. Furthermore, the manner in which Clark uses 
streamflow and/or rainfall data within a large region to establish flood fre- 
quencies for a small basin suggests a tendency to confuse the frequency of 
occurrence of an extreme event within the region with the probability of 
its occurrence on the small basin. 

In conclusion, next to streamflow measurements, rainfall data are logical 
the best index to flood frequencies. The idea of regional procedures is soun 
but there is no reason why such procedures based on the relatively short 
records of streamflow generally available should be reliable for the longer 
recurrence intervals. The authors’ procedure could be used to increase 
synthetically the length of streamflow records for individual basins to pro- 
vide a more reliable base for regional procedures that would yield much mo 


dependable results, especially for the longer recurrence intervals, than is 
now possible. 
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TURBULENCE CHARACTERISTICS OF THE HYDRAULIC JUMP® 


Closure by Hunter Rouse, T. T. Siao, 
and S. Nagaratnam 


| HUNTER ROUSE,! M. ASCE, T. T. SIAO,2 and S. NAGARATNAM?.—Com- 
ments upon this paper were of particular interest in that they reflected the 
‘views of investigators at five of the principal hydraulics laboratories now en- 
|gaged in research on turbulent flow. Significant points already considered 
during the preparation of the paper were emphasized, and others not con- 
sidered or at least not given their proper attention were very pertinently 
noted. As is always the case, the net result is a marked improvement over 
the original contribution. 

The writers recommend that everyone who has access to a glass-walled 
flume carry out the experiment described by Mr. Silberman - artifically 
varying the velocity distribution, if necessary, by the application of a coarse 
bed roughness. The formation of a standing eddy on a level floor is as con- 
vincing a demonstration of a basic boundary-layer principle as one could 
hope to find. Other hydraulic distortions of the jump proper have been dis- 
cussed in detail elsewhere, but they might well be mentioned here in this 
connection: those due to variation of channel cross section, particularly 
through the trapezoidal range; to variation of bed slope, as in a chute, bed 
curvature, as on a spillway bucket, or bed alignment, as at the end of an 
apron; and variation in wall alignment or curvature, as in flared or warped 
chute spillways. Each of these will have an influence upon the jump charac- 
teristics just as important in its way as that of bed roughness - but just as in- 
essential to the basic jump phenomenon. In conducting the investigation, the 
writers attempted to eliminate every possible extraneous factor - including 
initial velocity distribution - so that full attention could be given to the 
mechanism of turbulent expansion. Surely the introduction of a sonic throat 
and shock wave, as suggested by Mr. Silverman, would not have been an ef- 
fective step in this direction. 

Messrs. Peterka and Bradley introduced a welcome note into the literature 
in their immediate recognition of the usefulness of academic findings in 
clarifying the practical aspects of a related problem. The question of jump 
length has long been a matter of controversy, in part for lack of agreement on 
nomenclature and in part because the jump actually never ends but only 
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approaches asymptotically a state of depth, velocity, and turbulence equilibri- 
um. Designation of an end is thus purely arbitrary. Messrs. Peterka and 
Bradley have realistically adopted a criterion based upon bed velocity; 
Bakhmeteff and Matzke logically based theirs on depth; the writers con- 
veniently terminated their boundary curves where the effect of further curva- 
ture seemed to be inappreciable. Were residual turbulence (the academic 
choice) to determine the nominal end, its position would obviously depend 
upon the percentage arbitrarily chosen for the permissible turbulence re- 
mainder. Since this remainder continues to be appreciable in comparison 
with that of the subsequent uniform flow for some distance downstream, such 
a definition would produce by far the greatest nominal length. The discussors 
are right in their belief that length of roller and length of rapid expansion 
would approach equality as the Froude number became very large; if the 
residual turbulence were related to the overall loss, all definitions would 
then agree at the (unfortunately, purely academic) limit. 

Mr. Robertson, to select a minor comment from his discussion at this 
point, speaks of a normal or “proper” hydraulic jump as one which occurs in 
established open-channel flow. This recalls, of course, the discussion by Mr. 
Silberman. Messrs. Peterka and Bradley, on the other hand, may well think 
of a proper jump as one controlled by the baffle piers of a foreshortened 
apron. “Established” flow is actually no more definite than is the design of a 
pier and apron, for it can vary through wide limits with channel roughness 
and shape. For this reason Mr. Harleman’s detailed momentum analysis of 
the jump, and his application of the resulting equation to the three conditions 
described in the paper, represent supplementary material of basic import. 
The writers fear, however, that readers will ignore the equation and assume 
the plotted curve to indicate a general solution. The fact should hence be em- 
phasized that even in established flow there will be a different deviation from 
the elementary equation for every Reynolds number, relative roughness, and 
channel shape. 

To return to Mr. Robertson’s discussion, the writers can only report that 
there was no indication of bottom separation in either flume or duct for the 
condition of minimal boundary-layer growth that was sought experimentally. 
So far as the upper boundary effect was concerned, in the vicinity of the non- 
uniformity it was surely not one of shear; whereas a more refined method of 
measurement would undoubtedly have revealed a positive shear above the 
roller, as correctly expected by Mr. Robertson, comparison of its probably 
magnitude with that in the region of high velocity gradient below would not 
lead one to anticipate much change from the conditions shown. The effect of 
the boundary upon large-scale fluctuations, on the other hand, was admittedly 
considerable, and some indication of its magnitude is undoubtedly to be found 
in Mr. Hubbard’s measurements. 

Free-turbulence literature now distinguishes among three different re- 
gions of the turbulence spectrum: the fine-scale region, containing the eddie: 
in which most of the energy of the turbulence is concentrated; and the large- 
scale region of intermittency described by Mr. Robertson. In this sense the 
large eddies which are so obvious a part of the true hydraulic jump represen’ 
neither an appreciable concentration of energy nor an important source of 
dissipation. They can, however, play a considerable role in the process of 
mixing. At the juncture between the roller and the main stream the intermit. 
tent formation of these large eddies (often distinguishable in jump photo- 
graphs) is much like that observed at the edges of jets, wakes, and boundary 


2 DISCUSSION 101 


ayers. The free surface, which is the more noticeable zone of intermittent 
luctuation and that to which Mr. Robertson refers, has some of the same 
racteristics. However, as noted in the discussion, it also has two inter- 
related features which distinguish it: abrupt discontinuity of the liquid 
edium itself, and the stabilizing or restoring action of gravity. The air 
odel surely duplicated the essential freedom of motion in the lower zone of 
intermittency to what extent the fluctuations in this zone are actually con- 
trolled by the free-surface fluctuations remains to be determined. 
| Mr. Robertson correctly implied that the plotted dissipation functions were 
‘of a partially assumed rather than a completely measured nature, in view of 
ithe representation of the nine velocity derivatives by a single one through the 
‘relationships of isotropy. Whether this assumption requires only a high Rey- 
Bias number to be closely valid, or the components even then would deviate 
ito different degrees in different directions and different zones of flow, must 
jalso remain to future investigations for a conclusive answer. Surely the 
present measurements do not indicate the existence of much variation in 
validity from section to section, else the production-convection -dissipation 
sum would have departed more radically than it did from a zero value. It 
‘should also be remarked in this regard that such indications of small-scale 
‘isotropy as exist to date have been found under conditions of free turbulence, 
‘whereas deviations therefrom have all been noted in the vicinity of solid 
boundaries. 

Measurements of the type proposed implicitly by Mr. Robertson and ex- 
plicitly by Mr. Hubbard must eventually be made, although the writers shud- 
der to think of the painstaking labor involved in determining at many succes- 
sive sections - under conditions in which even the depth itself is very hard to 
define - what to date has been measured only for representative sections of 
flows which are either uniform or similar from section to section. Mr. Hub- 
bard apparently has at hand equipment that will accomplish some part of the 
measurements required - at least at a safe distance below the freely 
fluctuating surface - and he has demonstrated a knack of filling other needs 
as they arise. Still to be found are’a group of expert laboratory technicians 
with considerable time at their disposal, and an unusually patient and 
sympathetic sponsor! 
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SNOWMELT RUNOFF2 


Discussion by K. R. Wright 


: 
: 


_ KENNETH R. WRIGHT,! A.M. ASCE.—Snowmelt forecasting is usually 
nsidered as being limited to the mountainous regions of the country. In 
S paper the authors have shown that it can also be utilized in non- 
Ountainous areas. 

The method presented by the authors appears to be a sound approach to 
ae problem of snowmelt forecasting, and it does present a procedure for 
eveloping similar relationships for other non-mountainous areas. 

The writer has checked the reliability of the relationship 


Ye ATP NTS QUEL 


Where X= Me/é Potentia/ 
Y=.% of fPunoff atl 7: MP 


y applying it to three years of records in much the Same manner as a river- 
vel forecaster would use it. The three years chosen were 1953, 1954, and 
955, which were subsequent to the authors’ period of study. 

The quarter-day dry bulb temperature, dew point temperature, wind speed, 
nd sky cover for the Wausau, Wisconsin, station were used to calculate the 
larter-day melt potentials. The melt potentials were a summation of the 
Invection, condensation, and radiation melt potentials, as was stated by the 
ithors. 

Precipitation was also determined on a quarter-day basis throughout each 
‘the three spring seasons. The Medford and Marshfield, Wisconsin, sta- 
ons were used to determine the precipitation. These are the two stations 
3ed by the authors. The rainfall was added directly to the total quarter-day 
elt potentials. 

Daily discharges of the Big Eau Pleine River as given by the U. S. 
ecological Survey were corrected according to the “ice index theory,” given 
the paper, for the period of ice effect. The beginning of the melt period 
id the total inches of runoff were determined from these daily discharges. 

The meteorlogical data was then used as a river-level forecaster would 
e the data to determine future flows. The forecasted flows were then 
ecked against the actual daily river discharges. Due to the unusually small 
iow accumulation in 1954, however, it was necessary to resort to the 
thors’ Table VI to forecast the snowmelt hydrograph for 1954. 


Proc. Paper 1834, November, 1958, by J. H. Zoller and A. T. Lenz. 
Cons. Engr., Boulder, Colo. 
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The forecasted hydrographs are compared to the recorded hydrographs i 
Figures 1 and 2 for 1953, 1954, and 1955. In 1955 the forecasted and re- 
corded hydrographs are quite close for the first five days of the snowmelt, 
and then they diverge. The divergence is due to the fact that the snow on the 
ground had completely melted. On March 21 and 22, though, there was a nev 
snowfall which required ripening, and which then ran off on March 30 and 31 
as indicated by the hydrograph. The forecasted hydrographs were based on 
the runoff reaching the gaging station the day after the melt potential oc- 
curred. This is in accordance with the authors’ statement that there is a 
13-hour lag in runoff for the basin and also that most of the melt potential 
occurs in the second period of each day. 

The forecasted and recorded runoff for the three years was put into mass 
curves and also compared. These are shown in Figures 3 and 4. The hori- 
zontal portions of the forecasted mass curves are due to there being no mell 
potential for those specific days. This is illustrated for 1955 from March 1 
through March 19. 

In Figure 9 of the authors’ paper the relationship developed is presented 
graphically and compared with several years of reliable data. Figures 5 ant 
6 of this discussion give the relationship as determined for 1953 and 1955. 
It is shown that the equation which was developed fits well with the recorded 
runoff for these two years. 1954 was not compared because the very small 
amount of snow on the ground caused a distortion of the relationship. 

The summary of the melt conditions for the three years studied for this 
discussion is as follows: 


Length of Rain During Inches of Runoff 
Year Melt Period Melt Period Melt Period During Melt Peri 
1953 March 20 to 9 days 0.12 inches 4.17 inches 
March 28 
1954 March 24 to 2 days 0.58 inches 0.76 inches 
March 25 
1955 March 12 to 21 days 1.15 inches 3.77 inches 
April 1 


The Climatological Data bulletins report that in 1953 there were 15 to 20 
inches of snow on the ground before melting began. In 1954 there was only 
about one inch of snow, resulting in only a very small amount of runoff. Th 
1955 bulletin lists 3 to 4 inches of snow before melting, and then on March £ 
and 22 a snowfall raised it to 8 inches. 

From the writer’s investigation of the three years subsequent to the 
authors’ period of study, the following conclusions were drawn in regard to 
the Zoller-Lenz method of forecasting snowmelt runoff. 


1. The forecast of the daily discharges for 1953, 1954, and 1955 correlat 
satisfactorily with the recorded discharges on the Big Eau Pleine River. _ 


2. The relationship developed does not hold true for the runoff from un- 
usually small snow accumulations such as in 1954. 


3. With a small drainage basin the size of the Big Eau Pleine Basin, the 
runoff forecasts can only be made for one day in advance. The U. S. Weath« 
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reau’s five-day forecasts cannot alleviate this to any great extent because 
the amount of detailed meteorlogical data needed for the authors’ method. 


4. The determination of when the snow will attain ripeness and thus when 
itional melt potential will cause runoff is a problem which has not been 
lvyed by this paper. On small drainage basins such as the Big Eau Pleine 
Sin, the forecaster would probably miss the first day’s snowmelt runoff. 

» would, however, be alerted to the ripeness by this first day’s runoff and 
m begin his forecasts for the subsequent days. 


5. The “ice index theory” appears to be a sound approach to the problem 
eo aining discharges during the rising stage of the hydrograph. 


In summary, this paper on snowmelt runoff is a valuable contribution to 

published information on hydrology. The fact that this study was based on 

on-mountainous region indicates that the interest in snowmelt runoff is 
reading beyond the areas where it comprises nearly the entire supply of 
ater. 
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